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THERMAL AND ELECTRICAL PROPERTIES ARMCO 
IRON HIGH 


ABSTRACT 


Thermal and electrical conductivity, thermoelectric power vs. platinum, 
and thermal expansion Armco iron were determined the temperature range 
0°C 1000°C. All these properties show discontinuous change the 
transition iron, and change slope the Curie point. These measure- 
ments were carried out contribution co-operative determination 
thermal conductivity Armco iron high temperatures. 


INTRODUCTION 


The thermal and electrical properties iron high temperatures have 
been repeatedly measured the past. The agreement between the various 
experimental results the whole quite satisfactory. The one exception 
the thermal conductivity, particularly temperatures above the Curie 
point. far three experimenters have published results measurements 
thermal conductivity temperatures above 800° These are Maurer (1936), 
Powell (1939), and Lucks and Deem (1958). Powell made both comparative 
measurements for linear heat flow and absolute measurements for radial 
flow cylinder. Maurer made absolute measurements for linear heat flow 
and Lucks and Deem comparative measurements for linear heat flow. Oddly, 
and perhaps significantly, none the experimenters give estimates the 
error, but 1000° their results differ 27% the mean value. 

The situation much better below the Curie point. Three further deter- 
minations have been done the range all absolute and 
with linear heat flow, Powell (1934), Hattori (1937), and Armstrong and 
Dauphinee (1947). These authors give estimates their experimental 
error which are less than for the results Powell, 1.95% for Armstrong 
and Dauphinee, and 0.8% for Hattori. The maximum spread results, 
between Armstrong and Dauphinee and Hattori, the average 
30° and 9.7% 700° much larger than the sum the quoted errors. 
There is, course, the possibility that least part the discrepancy 
the experimental results, particularly lower temperatures, due the 
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variation the chemical composition the samples (Table I), but 
inconceivable that could account for the large spread high temperatures. 


TABLE 
Chemical composition the specimen Armco iron used experiments 


Percentage composition iron 


(by 
Powell 0.023 0.007 0.025 0.007 0.020 
Maurer 0.01 trace 0.011 0.040 0.05 
Hattori 0.015 0.010 0.032 0.003 0.030 
Armstrong and Dauphinee 0.039 0.067 0.024 99.807 


*As given by the distributor. 


resolve this conflict experimental results, Lucks the Battelle 
Memorial Institute proposed co-operative measurement thermal con- 
ductivity Armco iron. Four laboratories agreed participate: the National 
Physical Laboratory, Teddington, England, the National Bureau Standards, 
Washington, U.S.A., the Battelle Memorial Institute, Columbus, U.S.A., and 
the National Research Council, Ottawa, Canada. The Battelle Memorial 
Institute obtained the specimen Armco iron the shape in. diameter 
rods, annealed them for half hour 850° and distributed them the 
participants. This laboratory received long section marked No. 
and No. its ends. The sample for the experiment was cut out from the 
middle the section and X-rayed. The radiograph revealed cracks 
non-uniformities. The thermal and electrical conductivity and the thermo- 
electric power measurements were made this sample. The thermal expansion 
measurements were made another piece Armco iron. The method, 
results, and errors are discussed separately for each the properties the 
following section. 


METHOD, RESULTS, AND THEIR ERROR 


Thermal Conductivity 

Experimental Method 

When Battelle proposed the co-operative thermal conductivity measure- 
ment, they suggested aim the determination the conductivity the 
temperature range 1000° with experimental accuracy 3.5% 
better. However, the dimensions the sample supplied were not ideally 
suited such experiment. They necessarily forced the adoption 
linear heat flow method, without giving quite enough room make sure that 
the flow was rectilinear. The relatively small diameter the sample made 
difficult devise heat source that was truly uniform, particularly near 
the edges, and limited the temperature measurements the surface the 
sample, where the edge non-uniformities the heat source have the greatest 
When possible errors due these effects were estimated, was seen 
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that the conventional matched guard technique was not good enough 
ensure error less than The matching error is, course, 
dependent the geometrical ratios and thermal parameters the com- 
ponents employed the apparatus; for reasonable system designed around 
the specimen supplied, the error due single guard about 0.2% for perfect 
lateral guard and sample temperature matching three points, and increases 
about 0.5% for each degree lateral mismatch. One cannot expect 
match the sample and guard temperatures less than 1000° (see, 
for instance, the experimental results Powell (1934)), which puts the mis- 
match error alone about 2%. eliminate this error method was developed 
that does not depend temperature matching sample and guard. brief 
description the method here given with reference the actual experi- 
mental arrangement, shown Figs. and exact mathematical 
analysis given the Appendix. 

The sample, which two halves, and its environment shown Fig. 
There, (4) and (8) are the top and bottom parts the sample, each length 
7.156 and diameter 2.324 cm; (6) thin central heater between them, 
equipped with current and potential leads; (1) and (11) are the heat sinks. 
The contact surfaces between the sample, the central heater, and the heat 
sinks (e.g. (9) and (17)) were coated with high purity tin and sweated 
together. Each part the sample had three thermocouples (5) wedged into 
small holes (0.4 diameter and 0.5 deep) that were straight 
line 3-cm intervals, the central thermocouple being equidistant from the 
heater and the sink. The thermocouples, Pt/Pt-10% 0.25 diameter, 
were inserted into fine alumina tubes, taken halfway around the sample 
the level the point contact and then brought out through the bottom 
heat sink. 

Around the sample was alumina guard (7), with internal diameter 
4.96 cm. Five heaters (15) were wound the outside the guard, and 
six thermocouples (16) were positioned opposite the sample thermocouples. 
These thermocouples were taken halfway around the guard before they were 
brought out sideways through the radiation shields surrounding the guard. 
Grounded wires (14) were cemented the guard either side the thermo- 
couples minimize leakage from the heaters. Around the guard were several 
radiation shields supported between the end plates (3) and (13). End heaters 
(12) and (18) were wound the top and bottom sinks control the tem- 
perature the samples; these had several radiation shields around them ((2) 
and (10)). The space between the guard and sample was filled with fine 
alumina powder, and the whole arrangement was placed vacuum, shown 
Fig. Here (3) the water-cooled vacuum vessel, with the port (5) con- 
Hg. Both the top and bottom sinks had constriction them (2) limit the 
loss heat high temperatures, and had water cooling provided (4) 
withdraw heat lower temperatures. Pressure was applied vertically the 
sample spring (1) maintain good thermal contact temperatures 
above the melting point the tin the contact faces the sample. All the 
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Fic. Experimental arrangement the sample and guard. 
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Arrangement the sample guard system the vacuum vessel. 
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heater and thermocouple leads were brought out through the bottom the 
vacuum vessel through kovar seals. The heaters were controlled variacs 
powered from isolating transformers connected stabilized voltage source. 
The thermocouple e.m.f.’s were fed into 12-point recording potentiometer 
for rough measurements; thermal equilibrium the e.m.f.’s were switched 
over Wenner potentiometer and measured 1/10th microvolt. The 
vacuum, measured ionization gauge, was room 

the unmatched guard method the sample and guard are treated mathe- 
matically single system. The system must fulfil four conditions: 

(1) That the plane through the center the central heater and coplanar 
with plane zero longitudinal temperature gradient for the system. 

(2) That the medium between sample and guard have conductive mode 
heat transfer. 

(3) That the temperature each the sink ends the system uniform. 

(4) That the temperature variation the conductivities all the com- 
ponents the system small for such temperature differences may exist 
the system. 

Under these conditions the determination the thermal conductivity 
involves two sets temperature readings, with and without power developed 
the central heater. These readings were taken follows: The end heaters 
((12) and (18) Fig. and the guard heaters were adjusted until the sample 
was the desired temperature and the plane through the central heater was 
plane zero longitudinal temperature gradient. When thermal equilibrium 
was reached, the e.m.f.’s the six sample thermocouples and the six guard 
thermocouples were measured set order fixed time intervals. Two sets 
readings were taken, the readings for each thermocouple being minutes 
apart. (These fixed time intervals permitted the calculation correction 
for possible time variations temperature.) Then, the central heater was 
switched on, dissipating approximately watts. The power was measured 
with Weston Model 310 wattmeter, and was corrected for the power loss 
the wattmeter. The power the end heaters was reduced that the 
temperature the sink ends the sample did not change appreciably, and 
the end heaters were then adjusted that the central plane remained plane 
zero longitudinal temperature gradient. The heaters the guard remained 
untouched during this procedure. equilibrium, two sets temperature 
readings again were taken. Under these conditions the surface temperatures 
each half the sample are given 


where the distance measured along the axis the sample, corre- 
sponding the plane zero longitudinal temperature gradient, the 
length each half the sample, the difference between the tempera- 
tures measured the surface the sample height with the central heater 
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and off, and the gradient due the central heater set each half 
the sample, that 


(2) 


where the power input into the central heater, the cross-sectional 
area the sample, and the average conductivity the sample. The 
subscripts and refer the top and bottom half the sample. may 
well mentioned here that ranged from 10° C/cm, and that the 
average temperature the two parts the sample never differed more 
than 

shown the Appendix, the infinite series equation converges 
rapidly for K,/K, 100, where the sample conductivity and the 
conductivity the powder between the guard and sample. our experiment 
fore equation can written 


The temperature being measured three points, three unknowns could 
evaluated from equation These would normally and Ao. 
However, was kept very small, about 1°C 3°C, different 
calculating procedure was employed: was extrapolated from the 
measured the equation 


AT(z) 


This should give very good approximation the true the gradient 
the sample almost linear, and small, being the position 
the thermocouple nearest the sink. The calculated was then used 
determine and from equation and was subsequently calculated 
from modified form equation 


This calculating procedure giving well was employed not from 
any mistaken idea improving the accuracy but obtain estimate 
the correspondence the experiment with theory. The theory gives definite 
expressions for the terms the geometry the system and the con- 
ductivity ratios K,/K,. However, known only very approximately 
the agreement between theory and experiment can judged solely the 
ratios the and not their absolute value. 

can seen from the preceding description the method that the guard 
does not play important role the measurements. there merely 
help establish the plane zero temperature gradient across the system, and 
keep the temperature the whole system reasonably uniform, that 
the temperature variations the thermal conductivity the components 


the system small. 
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The actual experiment consisted four separate runs, made with three 
sets thermocouples and two different central heaters. The outside dimen- 
sions the two central heaters were the same, 2.2 diameter and 0.6 
thick. One them, however, was made out steel, very much like that 
Powell (1934). The other one was made out alumina, into the middle 
which was sintered fine coil Pt-20% wire going zigzags right 
across the heater. Two platinum disks were sintered onto each face the 
central heater, that could tinned the sample. Although the ceramic 
heater was the stronger the two mechanically, the temperature distributions 
across its faces were not good for the steel heater. 

The experiments were done the following order. First run: ceramic central 
heater, first set thermocouples. Second run: ceramic central heater, second 
set thermocouples, position the two parts the sample interchanged. 
Third run: steel central heater, second set thermocouples, sample positioned 
second run. Fourth run: steel central heater, third set thermocouples, 
sample position second run. The sample, sink, and central heater surfaces 
were retinned before each run. The thermocouple separations and their actual 
distances from the sink were measured before and after each run. The thermo- 
couple separations were measured passing d-c. current through the 
sample and comparing the e.m.f.’s across the thermocouples those across 
knife-edge standard. The standard was placed against that half the 
sample for which the thermocouple separation measurements were made 
three different azimuthal positions relative the thermocouples. The separa- 
tion the knife edges the standard was approximately its exact 
value being measured with travelling microscope before and after each set 
measurements. The distance from thermocouple sink was also measured 
electrically, the d-c. e.m.f. between any one the thermocouples and 
needle probe placed the joint sample and sink being compared with the 
e.m.f. across the thermocouples. From the known dimensions the sample 
and central heater the distance the thermocouples from the plane zero 
longitudinal temperature gradient could calculated. All the electrical 
measurements were done with the complete system assembled, except that 
the ceramic central heater was used, copper plug would substi- 
tuted for it. The sample dimensions were also checked before and after every 
run. 

Experimental Results 

The results obtained the four runs are shown Fig. together with 
those Powell, Maurer, and Lucks and Deem. The results shown have 
had two corrections applied them. The first was correction for the non- 
uniform power distributions the central heaters. This was done measuring 
the surface temperature variation each heater along two mutually per- 
pendicular diameters, under conditions approximating those the experiment. 
These measured temperature distributions, 6), were fitted Bessel 
functions 


a 
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where measured with respect the position the sample thermocouples, 
constant, and the m’s are chosen that being the 
radius the central heater. Usually four terms with were sufficient 
obtain good fit. Under the assumption heat loss from the surface 
the sample, the sample surface temperatures will given 


sinh 


from which the required correction could calculated. These were: 


ceramic heater, top surface 1.027G 
ceramic heater, bottom surface 1.007G 
steel heater, top surface 1.001G 
steel heater, bottom surface 


where the temperature gradient calculated from equation and 
the corrected gradient. The correction factors are unfortunately not very 
accurate, due the difficulty obtaining meaningful surface temperature 
distributions. Successive traverses along one line differed magnitude 
equivalent 10% GL, which, the correction factor represented 

The other correction deals with deviations from thermal equilibrium. 
the experiment, the thermocouple e.m.f.’s were recorded, that their time 
variations were easily observed. was noticed that after every change the 
heaters the thermocouples would quickly settle down relative each other, 
usually the order minutes, while frequently there was slow, 
uniform drift all the thermocouples which, high temperatures, was 
difficult eliminate, but very easy correct for. The correction calculated 
for linear heat flow with the suitable boundary conditions, 


where the corrected and the uncorrected gradient, the average 
rate temperature change for each half the sample, the subscripts and 
referring the readings with the central heater and off, respectively. 
the diffusivity iron, and the time interval between the readings 
fixed time intervals, the correction was easy evaluate. The largest correction 
made was 0.7%, and the average 0.06%. 

can seen Fig. that the corrected experimental points lie two 
distinct curves, one for each the central heaters, while there significant 
difference associated with the different sets thermocouples used. The results 
for the ceramic and steel central heaters overlap the range 400° 
800° where the difference between the two sets (3.4+0.3)% their 
average. Equal weight was given the two sets results, for although there 
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were more points taken with the ceramic central heater, the correction that 
had applied them due the non-uniformity the power input was 
larger than that applied the points taken with the steel heater. (As matter 
fact, the two sets results were separated only 1.4% when uncorrected, 
the set for the steel central heater being higher than that for the ceramic 
central heater. The power distribution correction pushed the ceramic heater 
results further down approximately 1.7%, and the steel heater results 
further approximately 0.4%.) The average value for the thermal con- 
ductivity taken equal the mean the two sets the temperature 
range 400° 1.017 times the value for the ceramic central 
heater the temperature range 400° and 0.983 times the 
value for the steel central heater the range 800° 1000° The average 
value computed shown Fig. and given Table together with 


200 400 600 800 1000 1200 


Fic. Thermal conductivity Armco iron function temperature: experi- 
mental points the first run; experimental points the second run; experimental 
points the third run; experimental points the fourth run; average value for this 


the values the other experimenters. The values listed Table are not 
corrected for the thermal expansion iron. 


Experimental Error 

The uncertainties the determination the thermal conductivity will 
discussed with reference equations and There are two uncertainties 
associated with the determination the thermal conductivity, asa function 
temperature: the probable error the determination and the probable 
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TABLE 
The thermal conductivity Armco iron given various observers 


Thermal conductivity, watts/°C 


Lucks and 


Powell Maurer Hattori Dauphinee Deem Laubitz 


*Taken from graph. 


error the measurement the temperature which supposed 
correspond. this experiment, the average the readings the four extreme 
thermocouples the two halves the sample, measured with the central 
heater on, was taken the temperature that particular point. Thus, 
individual errors the thermocouples tend average out, and the main error 
arises from the composition that particular lot wire from which the 
thermocouples were made. Deviations due wire composition are usually less 
than 5°C and and therefore the uncertainty the temperature 
should not exceed 0.5%. Effects due contamination, cold work, etc., the 
thermocouples appear negligible view the agreement the results for the 
various sets thermocouples, and for points taken the same temperature 
the beginning and the end run. 

The uncertainty arises from uncertainties the measurement 
and the power input, was measured with calibrated precision watt- 
meter, with error less than 0.2% full scale. For the power input used 
the experiment the wattmeter read the midpoint its scale, and therefore 
the uncertainty less than 0.4%. the cross-sectional area, was measured 
before and after each run. The maximum change observed during run 
was 0.26%, which introduces uncertainty not more than 0.13%. 

The uncertainty the measured temperature gradient, more difficult 
evaluate. arises due errors and uncertainties the parameters entering 
into equation directly, and from the experimental non-conformity the 
conditions assumed deriving equation Into the first category fall the 
uncertainties the temperature differences and the thermocouple 
separation The uncertainties the absolute values are not 
important, the correction term and are less than 
The being small temperature difference given thermocouple, 
should have fractional error larger than the fractional error the tem- 
perature measurement, that 0.5%. small compared 
the uncertainty only slightly larger than the uncertainty 
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about 0.6%. The thermocouple separation was measured with 
precision 0.2%, but was found that the numerical value would change 
during run, the average, 0.8%, presumably due changes the 
tinned contacts the sample. The uncertainty therefore the 
order 0.6%. 

There are several factors contributing the deviation the experiment 
from theory. First, there the finite size the central heater. This, however, 
turns out negligible; according simplified calculation, the ceramic 
central heater affects the temperature distribution the sample less than 
0.01%. Then there the non-uniformity heat input the central heater. 
correction has been made for this, but the correction factor itself accurate 
only 0.5%. There also the non-uniformity the contacts the sample 
with the central heater and the sinks. has been found that with all the 
interfaces freshly tinned and carefully joined, the measured thermocouple 
separation function the azimuthal position the standard, the variation 
being the order 1.2%, and that the change such azimuthal variation 
during run the same magnitude. taking the thermocouple separa- 
tion the average value for several azimuthal positions the standard, the 
effect the variation wiped out. The change the variation, however, 
does introduce uncertainty half the observed change, 0.6%. There may 
be, course, also radial variation the quality the contacts; the magni- 
tude the uncertainty that this would introduce not known, but one may 
expect comparable the uncertainty due the azimuthal variation, 
that approximately 1.2%. 

Finally, there will deviation experiment from theory due such factors 
non-uniformity thermal conductivity the powder separating guard and 
sample, non-uniformity sink temperatures, lack radial symmetry the 
guard, etc. These are the most difficult calculate, but estimate can 
made from the ratio calculated from equation Figure gives the 
experimental value where 


A,/GL 


while, according theory, all are positive, and for the 
above value should the order can therefore say that 
the uncertainty least the order experimental, and probably 
more, the order twice which amounts course, part that 
uncertainty probably caused factors considered earlier, which does 
complicate the assessment the total probable error. However, was felt 
that the square root the sum squares the individual uncertainties would 
probably fair indication the accuracy the measurement. Calculated 
that way, the over-all uncertainty 2.1%. The 0.5% uncertainty the 
temperature associated with can interpreted further error 
increasing the latter about 2.2% 100° and 2.5% 800° The observed 
difference 3.4% between the results for the two central heaters therefore 
well within the estimated uncertainty 
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200 400 600 800 1000 


The over-all experimental uncertainty should not confused with the 
precision the results for each run. Thus, for the ceramic central heater, the 
maximum deviation point from the smooth curve drawn through all the 
points was 1.9%, and the r.m.s. deviation from that curve only 0.48%. For 
the steel heater, the maximum deviation from the curve was 0.9%, and the 
r.m.s. deviation 0.47%. 


2.2. Electrical Resistivity 

The electrical resistivity was measured during and after the last run with the 
steel central heater. Measurements were done with the sample thermal 
equilibrium, and the heaters adjusted that the temperature throughout the 
sample did not vary more than The measurement involved passing 
d-c. current through the sample, and comparing the e.m.f. across the two 
platinum legs the outside thermocouples each part the sample with 
the e.m.f. across standard resistor series with the sample. Six measure- 
ments were done each temperature for each half the sample, the current 
being reversed after every measurement. The average value any tempera- 
ture for both halves the sample shown Fig. The r.m.s. deviation from 
the curve drawn through these points 0.59% for one half the 
sample, and 0.67% for the other, while the r.m.s. difference between the 
points the two halves the sample 0.69. This indicates that there 
consistent difference between the results for the two halves. 

The average value the resistivity given Table III. Also given there 
are Powell’s results for his sample and, for purposes comparison, the results 
obtained high purity iron wire Wallace, Sidles, and Danielson (1960). 
The average value the resistivity was used the calculation the Lorentz 
function, which shown Fig. together with the Lorentz function 
calculated Powell for his specimen. 

The uncertainties the values assigned for the electrical conductivity come 
from the error the determination the ratio the e.m.f. measured across 
the sample the e.m.f. the standard, and from the uncertainties the 
cross-sectional area the sample, the separation the thermocouples, and 
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200 400 600 800 1000 


Fic. The electrical resistivity and Lorentz function Armco iron: experimental 
points the electrical resistivity, Lorentz function calculated from the results this 
experiment; Lorentz function obtained Powell. 


TABLE III 
Electrical resistivity and the Lorentz function Armco iron 


Electrical resistivity, Lorentz function (V/°K)? 
Armco iron High purity 
Laubitz Powell iron Laubitz Powell 
100 15.6 15.1 2.7 2.76 
200 22.9 21.89 2.83 2.93 
300 32.0 31.4 30.97 2.95 3.01 
400 42.9 43.1 41.82 3.02 3.06 
500 55.2 54.39 3.08 3.08 
600 70.3 69.9 3.09 3.08 
700 87.9 87.0 85.85 3.07 3.06 
800 106.0 105.5 104.33 3.01 2.93 
900 112.0 110.78 
1000 115.3 115.4* 114.49 2.74 


*From graph. 


the temperature determination. Each e.m.f. was measured with accuracy 
0.6%, that the average the ratio should accurate about 
0.5%. The other three uncertainties are 0.2%, 0.6%, and 0.5% respectively, 
giving over-all uncertainty about 1.5% for the results each half the 
sample, and about 1.0% for their average. 
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The uncertainty the Lorentz function the combination the un- 
certainty the resistivity and thermal conductivity, about 3%. 


2.3. Thermoelectric Power Armco Iron vs. Platinum 

The e.m.f. between the platinum wires the three thermocouples each 
half the sample were determined simultaneously with every measurement 
temperature during the thermal conductivity experiment. These differential 
measurements were used checks the internal consistency the 
Pt/Pt-10% thermocouples. Furthermore, subtracting the 
taken with the central heater and off, and dividing the corresponding 
difference the temperature differentials the thermoelectric power 
iron vs. platinum was calculated with accuracy 0.5%. The results 
obtained for the two outside thermocouples each half the sample are 
shown Fig. and their average value read off the curve given 
Table IV. The interesting feature these results the change slope the 


20.0 


. 
. 
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200 400 600 800 1000 


Fic. The thermoelectric power vs. platinum and the thermal expansion Armco iron: 
experimental points for the thermoelectric power vs. platinum, the two halves 
the sample; experimental points the total thermal expansion, AL/L; total thermal 
expansion calculated from data given the Metals Handbook. 


Curie point, and the 30% drop the transition. There systematic 
difference between the results for the two halves the sample, which good 
indication the uniformity the iron supplied for the test. The r.m.s. devia- 
tion all the points from the curve drawn through them 0.4%. 
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TABLE 


Thermoelectric power vs. platinum and thermal expansion 
Armco iron 


Thermal expansion AL/L 
Thermoelectric 
vs. platinum Laubitz Metals Handbook 


200 14.10 25.4 
300 10.82 39.9 
400 8.98 55.4 
500 9.12 71.5 
600 11.34 87.6 
700 15.32 104 
800 19.34 121 
850 19.83 128 
940 14.24 112 

1000 15.20 125 
1100 147 


2.4. Thermal Expansion Armco Iron 

The thermal expansion measurements were made small piece Armco 
iron, the expansion being measured relative sintered alumina. The purpose 
the measurement was find the relative displacement the sample thermo- 
couples and the guard thermocouples the thermal conductivity experiment. 
was found, however, that changing measurements from relative 
absolute, using the values Ebert and Tingwaldt (1936) for the expansion 
alumina, very good agreement obtained with values given the Metals 
Handbook (1948) for iron the temperature range 700° was 
therefore decided include the results the expansion measurements this 
report. 

Figure shows the experimental results for the total thermal expansion, 
AL/L, and Table gives the average values read off the curve. could 
expected, there sudden drop the expansion the transition, 
but the observed decrease slightly larger than calculated from crystal- 
lographic data. The observed change lattice parameters implies density 
change 0.4%, while the change expansion observed here corresponds 
density change 0.7%. 

The accuracy the expansion measurements low; judging from the scatter 
the experimental points, which represent four independent runs, the probable 
error the expansion, about 10%. 

The values given Table can used correct the thermal conductivity 
and electrical resistivity for the effect thermal expansion iron. The 
correction factor simply for the electrical resistivity, and 
(1—AL/L) for the thermal conductivity. 


DISCUSSION RESULTS 


The results the thermal conductivity measurements reported here are 
very good agreement with those Powell, the range from room temperature 
the Curie point. This close agreement shown not the thermal conduc- 
tivity itself, but the correspondence the Lorentz functions (Fig. 5). One 
expects that differences the iron the samples will affect the Lorentz 
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function much less than either the thermal electrical conductivity, long 
the electrons that are mainly responsible for the heat transfer. The 
comparison the Lorentz function implies therefore that about 3/5 the 
discrepancy the thermal conductivity room temperature due 
basic differences the constitution the sample, presumably the 0.09% 
difference their purity. The remaining discrepancy, room temperature, 
presumably due experimental error, although part may due toa 
further effect the impurities their scattering the phonons. can 
expected, the effect the impurities decreases with temperature, and the 
two thermal conductivity curves gradually approach each other, meeting 
the Curie point. 

Comparison with other results not satisfactory. There little that can 
said about the results Hattori except that his experimental error 
probably much larger than his estimate it. Powell (1936) has already shown 
that Hattori’s method temperature measurement seriously affected the heat 
flow, which any case was probably far from the rectilinear flow assumed 
the calculations. Maurer’s results were obtained using method open the 
same criticism Hattori’s, and the good agreement with Powell’s results 
lower temperatures probably fortuitous. The only serious difficulty 
presented the results Armstrong and Dauphinee. They have performed 
what appears very good experiment, and yet they obtained the highest 
conductivity for sample least purity. The difference between their and 
our results room temperature, where the combined estimated error 
only 4%. 

temperatures above the Curie point the results the various observers 
spread much that impossible arrive any correlation. Our results 
differ from previous ones not only absolute value, but also the shape the 
temperature dependence curve. This seen enlarged scale Fig. 
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Fic. conductivity and electrical resistivity Armco iron the vicinity 
the Curie point (768° and the transition (900° C-910° C). Experimental points 
for the electrical resistivity and thermal conductivity respectively, uncorrected for thermal 
expansion; the same points corrected for the thermal expansion. 
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where the thermal conductivity and the electrical resistivity are shown after 
correction for the thermal expansion iron. seen that the experimental 
points for the electrical resistivity definitely imply discontinuous jump 
about 1.8% the transition. Such jump not surprising considering that 
other properties, such thermoelectric power, thermal expansion, and 
specific heat, also show discontinuities the transition. also not incon- 
sistent with Powell’s results, which display peculiar little wiggle that 
temperature region. Other observers have also reported unusual behavior 
the electrical resistivity the transition (see, for instance, Wallace, 
Sidles, and Danielson (1960)), but nothing definite our results. Any 
discontinuity the electrical resistivity should have corresponding dis- 
continuity the thermal conductivity, that temperature region about 
85% the heat carried electrons. fact, shown Fig. discon- 
tinuous curve can drawn through the experimental points for the thermal 
conductivity with slightly better fit than obtained with continuous curve. 
The discontinuity about 1.7% the right order magnitude. 


SUMMARY 

Thermal conductivity, electrical resistivity, and thermoelectric power vs. 
platinum were measured for single piece Armco iron the temperature 
range 1000° The thermal conductivity was measured new 
unmatched guard method with estimated uncertainty 2.5%, and 
very good agreement 800° with results obtained Powell. All the 
properties, including thermal expansion measured different piece Armco 
iron, indicate greater lesser discontinuous change the transition 
iron, and change slope the Curie point. 


APPENDIX 
Mathematical Description the Unmatched Guard Method Measuring 
Thermal Conductivity 


Figure represents schematically each half the sample, together with the 
powder between the guard and sample. and are the thermal conduc- 
tivities the sample and powder, assumed constant and uniform, the 
radius the sample, the inner radius the guard, and the length each 
half the sample. The assumed boundary conditions are: 


(1) radial symmetry, 


The last condition imposes restriction f(z). 
Under these conditions, the solution Laplace’s equation 
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r=b r=a 


Fic. representation half the guard-sample system. 


where 
4 To (ant )K (and) + 11 (na) K 0 (arn) 
To (an@ )K (ana) + pli (ana o(a,a) 
and 


L 


and are the associated Bessel functions. 

assume that wafer-thin central heater located and extending 
=a, which, when switched on, establishes the sample uniform tempera- 
ture gradient (—G) for (due the symmetry, that gradient zero for 
0). further assume that the switching the central heater does not 
affect f(z), but may change 7(L) define 


where the superscripts and refer temperatures with the central heater 
and off, then 
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where 


Figure shows the functions for the geometry used 
the experiment. has been measured air, where approximately 


Fic. Theoretical values the parameters function for the geometry 
used the experiment. 


mw/°C cm. vacuum expected less than that, and therefore 
the infinite series equation 

course, the true temperature gradient established the sample 
that fraction the heat input which flows along that sample. only the total 
input into the central heater, measured, then 
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where za? and the subscripts and refer the top and bottom half 
the sample. 

may mentioned here that condition (4) unnecessarily 
when the central heater switched on. Equation A.1 then becomes 


where for one half the sample and for the other 
half. Equation A.2 therefore remains unchanged 
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THE MAGNETIC SUSCEPTIBILITIES SOME DIAMAGNETIC 
ALLOYS 


THE PRIMARY SOLID SOLUTIONS CADMIUM, INDIUM, TIN, AND 
ANTIMONY COPPER AND SILVER; ZINC, GALLIUM, 
GERMANIUM, AND ARSENIC SILVER; COPPER 
SILVER AND SILVER COPPER! 


ABSTRACT 


The results the measurements the variation with concentration the 
magnetic susceptibility the solid solutions the following systems are pre- 
sented: copper with silver, cadmium, indium, tin, and antimony solutes; 
silver with copper, zinc, gallium, germanium, arsenic, cadmium, indium, tin, 
and antimony solutes. The density states for silver calculated using 
the approximation nearly free electrons. found that the density states 
decreasing electron atom ratio 1.0. further found that the 
method used earlier paper discuss the magnetic properties the a-solid 
solutions formed with copper zinc, gallium, germanium, and arsenic 
applicable the systems the present paper. 
concluded from the magnetic measurements that the Fermi surface touches 
the Brillouin zone face both silver and copper and that the rate 
decrease the density states with electron atom ratio given approxi- 
mately the approximation nearly free electrons. The solid solutions the 
system are anomalous. 


INTRODUCTION 


The results the measurements the variation the diamagnetic sus- 
ceptibilities with concentration the solid solutions formed with copper 
solvent the four metals following the periodic table were presented 
earlier paper (Henry and Rogers This paper will referred 
hereinafter was found that the rate which the alloys become dia- 
magnetic increases with the valency the solute. The increase diamagnet- 
ism was ascribed two major sources: (1) the diamagnetism the solute 
ion and its localized electrons, where the difference valency between 
the solute and the solvent; (2) the decrease the net paramagnetic suscepti- 
bility the conduction electrons which results from decrease the density 
states the Fermi surface with increasing electron density (Jones 1937). 

was thought desirable investigate: (1) the influence the same series 
solutes from the [Vth period another monovalent solvent, determine 
whether not the change susceptibility did depend part upon the solvent; 
(2) the influence the solutes from the Vth period with the corresponding 
valency, determine whether not the apparent valency effects existed 
there well. 
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The measurements appear indicate that the contribution the polyvalent 
solute ion the magnetic susceptibility the alloy fairly independent 
the solvent and that there additional contribution which proportional 
the electron atom ratio, with the constant proportionality depending 
upon the solvent and being independent the solute. 

order make quantitative calculation the observed initial slope 
the magnetic susceptibility, concentration curves, the method suggested 
was necessary calculate the density states for silver the method 
suggested Jones (1937). The results this calculation are presented 
Section 


PREPARATION ALLOYS 


The approximate purity, iron content atomic parts per million, and the 
source supply the metals used the preparation the alloys are: 

Copper, greater than 99.999%, 0.2 p.p.m., purchased from the American 
Smelting and Refining Co. 

Zinc, 99.999%, less than 1.0 p.p.m., donated the Consolidated Mining 
and Smelting Co. 

Gallium, 99.99%, 2.0 p.p.m., purchased from the Eagle-Picher Co. 

Germanium, 99.99%, 1.0 p.p.m., purchased from the Eagle-Picher Co. 

Arsenic, 99.95%, less than 1.0 p.p.m. iron, 0.05 wt. antimony, distilled 

Cadmium, 99.99%, less than 1.0 p.p.m., donated the Consolidated 
Mining and Smelting Co. 

Indium, 99.99%, 1.0 p.p.m., donated the Consolidated Mining and 
Smelting Co. 

Tin, 99.99%, 1.5 p.p.m., purchased from Johnson Matthey and Co. Ltd. 

Antimony, 99.99%, 5.0 p.p.m., purchased from Mackay, Inc. 

10. Silver, 1.0 p.p.m., electrolytically refined (Henry and Rogers 
1956a). 

The latter paper will referred hereinafter The iron analysis only 
has been made this laboratory. 

Although some the alloys were prepared described paper about 
half were melted and chill-cast apparatus which has been described 
elsewhere (Craw and Henry 1956). The alloys prepared this apparatus were 
superior the earlier ones that both lateral and vertical segregation were 
considerably reduced. Quartz crucibles were used all cases for melting 
this latter apparatus except for those alloys which contained gallium and 
indium; these cases either purified graphite alumina crucibles were 
employed. The alloys were cast atmosphere purified argon. After 
heat treatment, specimen 7.5X20 was machined from the ingot. The 
specimens were acid-cleaned before measuring remove iron contamination 

The alloys were heat-treated remove local inhomogeneities the con- 
centration and ensure single phase, face-centered cubic, solid solution. 
The details the heat treatment were varied according the requirements 
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the alloy system. brief description the relevant portions the phase 
diagrams and the heat treatments follows. 


For these four systems the boundaries determined Owen 
and Roberts (1939) and Hume-Rothery, Mabbott, and Evans (1934) are 
good agreement. The limits the a-solid solutions are: 
39.9 at.% 600° 42.6 at.% 300° 19.5 at.% 600° 
19.1 at.% 300° 10.9 at.% 9.8 at.% 

The alloys these systems were annealed 550° purified argon 
atmosphere pyrex containers, for days and then air-cooled prior 
machining. 


The limits the a-solid solutions these systems are: 
(Owen and Pickup 1933), 1.7 at.% 549° 0.28 at.% 300° copper- 
indium (Owen and Roberts 1953), 10.9 at.% 575° 6.8 at.% 470° 
(Raynor 1944), 15.8 at.% 550° 1.2 at.% 200° 
antimony (Hume-Rothery, Mabbott, and Evans 1934), 5.6 at.% 600° 
2.8 at.% 200° 

The alloys these systems were heat-treated prior machining follows: 
copper-silver, weeks then quenched iced water followed 
hours 850° and then quenched iced salt water; 
days 550° and then quenched rapid quenching furnace special 
design (Owen 1944); weeks 550° and then quenched 
iced water; weeks 550° and then quenched iced water; 
weeks 575° and quenched the rapid quenching 
furnace. all cases the annealing was carried out pyrex containers 
purified argon atmosphere except for the high temperature annealing and 
quenching the alloys for which quartz containers were used. 


The limits the a-solid solutions the systems are: silver-copper (Chaston 
1953), 14.0 at.% 779° 0.42 at.% 200° silver-zinc (Owen and 
Edmunds 1938), 33.6 at.% 39.7 at.% 331° 
(Owen and Rowlands 1940), 17.8 at.% 615°C, 14.9 at.% 300° 
which are lower than the results Hume-Rothery, Mabbott, and Evans 
(1934); silver-germanium (Owen and Rowlands 1940), 7.92 at.% 575° 
1.5 at.% 270° (Owen and Rowlands 1940), 8.5; at.% 

The alloys these systems were heat-treated prior machining follows: 
silver-copper, weeks 600° then quenched iced water followed 
hours 800°C and then quenched iced salt water; 
days 600° and then days 600° and then 
air-cooled; silver-germanium, days and then quenched the 
rapid quenching furnace; days 600° and then quenched 
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the rapid quenching furnace. The silver-copper alloy measured 717° 
was given the same heat treatment the ones measured room tempera- 
ture. addition was maintained 717° for hours prior measurement. 
change occurred the susceptibility during the measurement. 

The alloys highest concentration the following systems were examined 
micrographically for the presence phase: 
silver—arsenic, and phase was found. 

The analysis the silver-cadmium, and alloys 
was carried out measuring the difference between the lattice parameter 
pure silver and that the alloy and making comparison with the results 
Owen and Roberts (1939). The same method was applied the silver- 
zinc system using the results Owen and Edmunds (1938). All other alloys 
were chemically analyzed. all cases sample was taken from near 
possible both the top and the bottom the specimen. The error recorded 
with the analysis the probable error the mean value the composition 
determined from the two ends the specimen. measure homogeneity 
rather than indication the error the analysis. The estimated error 
all the analyses +0.05%. 


SUSCEPTIBILITY MEASUREMENT 


The apparatus, the method measurement, the errors susceptibility, and 
the method correcting for ferromagnetic iron were discussed earlier 
papers and All the measurements were made room temperature, 
23+2° except for two 717° pure silver specimen and 
alloy containing at.%. The apparatus used for these measurements 
will described later paper. 

One atomic part per million iron solution copper changes the suscepti- 
bility about 0.1% (paper view the purity the copper, the error 
from this source the order the experimental error and correction was 
applied the copper alloys. The solubility iron solid silver reported 
approximately 5.0 p.p.m. (Tammann and 1930). was not 
possible therefore determine the influence iron solution silver 
the susceptibility silver. assumed that the same order 
magnitude copper then, view the purity the materials, the error 
from this source may neglected. Another difficulty arose, however, 
connection with the silver alloys, presumably associated with the low solubility 
iron silver. There was tendency for the silver alloys have larger 
differences between and x,, where and x,, are the susceptibility 
maximum field and the extrapolated susceptibility respectively. Some points 
for which this difference was greater than 5.0% fell off the curve, outside 
the limits experimental error, and were discarded. There were points, 
however, for which the difference was less than 5.0% which did not fall 
the curves. Only those points for which the difference was 5.0% less have 
been reported. 

The results the measurements for the solid solutions are shown 
Tables II, and III, and Figs. and 
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TABLE 
Magnetic susceptibility solid solutions; silver with solutes from the Vth period 
Mean Mean Mean 
Alloy composition (c.g.s. (c.g.s. (c.g.s. 
(silver) (at.%) per gram) per gram) per gram atom) 
system 
1.17+0.07 —0.1809 —0.1825 —19.70 
1.24+0.01 —0.1787 —0.1812 —19.56 
3.88+0.04 —0.1824 —0.1838 —19.86 
4.09+0.01 —0.1825 —0.1840 —19.89 
6.34+0.09 —0.1854 —0.1857 20.09 
12.60+0.09 —0.1920 —0.1924 
19.20+0.07 —0.2012 —0.2020 —21.89 
26.26+0.11 —0.2074 —0.2115 
30.40+0.08 —0.2236 —0.2239 24.47 
system 
1.00+0.02 —0.1834 —19.80 
3.08+0.07 —0.1893 —0.1906 —20.60 
6.04+0.04 —0.1962 —0.1966 —21.29 
8.62+0.14 —0.2042 —22.25 
10.56+0.12 —0.2085 —22.71 
14.01+0.17 —0.2190 —0.2191 
17.06+0.04 —0.2272 24.88 
system 
—0.1868 —0.1874 20.24 
2.42+0.01 —0.1940 —0.1954 —21.14 
5.82+0.07 —0.2078 —0.2085 22.62 
7.30+0.11 —0.2139 23.39 
system 
1.05+0.01 —0.1853 —20.38 


The initial slopes the susceptibility, concentration curves were deter- 
mined fitting either straight line parabola, the method averages, 
the data (Table IV). 

The increase the rate change the susceptibility with increasing 
valency apparent Figs. and the three series Ag(V), Cu(V), and 
Ag(IV). The Roman numerals parentheses indicate the period the solute. 
this respect the behavior these three series the same that found for 
the series paper more detailed discussion the results 
given Section 

Other measurements the magnetic susceptibility the same solid solutions 
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TABLE 
Magnetic susceptibility solid solutions; copper with solutes from the Vth period 
Mean Mean Mean 
Alloy (c.g.s. (c.g.s. (c.g.s. 
(copper) (at.%) per gram) per gram) per gram atom) 
Copper-silver system 
1.16+0.10 —0.0872 —0.0878 
1.49+0.00 —0.0879 —0.0880 —5.65 
system 
Copper-indium system 
—0.1014 —0.1018 —6.66 
Copper-tin system 
2.60+0.02 —0.1012 —6.58 
system 
0.90+0.00 —0.0929 —0.0929 
1.84+0.07 —0.0997 —6.45 


are follows. Childs and Penfold (1957) measured one alloy the 
cadmium system containing 30.5 at.% cadmium. correction for ferromag- 
netic impurity was made extrapolating against 1/H. Their value for pure 
silver and for the one alloy are included Fig. Stephens and Evans (1936) 
measured one alloy the solid solution region the system 
containing 2.7, at.% antimony. not stated that these measurements were 
corrected for ferromagnetic impurity. Their value for pure silver and for the 
one alloy are included Fig. Endo (1925) measured one alloy the solid 
solution region the system containing at.% tin. not 
stated that correction was applied for ferromagnetic impurities. The value 
for pure copper and for the one alloy are included Fig. all cases the 
slopes agree roughly with those found the present work. The differences 
the absolute values are presumably due differences calibration and 
some cases differences the purity materials. Meyer and Weiner (1957) 
measured the temperature dependence several alloys the solid solution 
region the silver-zinc system. Comparison was not meant made 
among the alloys and thus their results have not been included. 
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TABLE III 
Magnetic susceptibility solid solutions; silver with solutes from the period 


Mean Mean Mean 
(silver) per gram) per gram) per gram atom) 


system 


3.02+0.00 —0.177 —0.1789 —19. 


717° 
—19.32 
system 
2.10+0.05 


—0.1750 


Silver-gallium system 


.46+0.00 —0.1873 
—0.1922 
—0.192 


Silver-germanium system 


1846 
.78+0.00 
1864 
2.38+0.02 —0.192 
—0.2023 


system 


.49+0.00 —0.1923 0.1937 
—0.2039 
5.26+0.01 —0.2126 


TABLE 


Calculated and observed slopes 


Experimental 


Slope 
difference difference 
Solute per gram atom per mole 


—16 

—18. 

26. 

35. 

—30.6 (2) 


*Numbers in parentheses obtained by using ZX 14 X 10-6 instead of ZX 17X 10-6. 


1831 
—19.75 
—20.40 
20.65 
20.80 
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Fic. susceptibility solid solutions; silver with solutes from the period. 


DENSITY STATES FOR SILVER 


Jones (1937) calculated the density states for copper method based 
the approximation nearly free electrons. result the hyperbolic 
form the surfaces constant energy the region the (III) zone face, 
the peak the density states function occurs electron atom ratio, 
e/a, approximately 1.0. Thus, when the electron atom ratio increases, 
assumed upon the addition polyvalent solute, the density 
states should decrease. attempting explain the magnetic properties 
the solid solutions the period solutes copper (paper the possi- 
bility decreasing density states was found convenient. 

Slater (1936) calculated the density states for copper using the data 
Krutter (1935) which was obtained using the cellular method. electron 
atom ratio 1.0 the density states still rising. Howarth (1953, 1955) 
when comparing the results the calculations the eigenvalues and eigen- 
functions copper calculated the cellular method and the method 
augmented plane waves, emphasizes the sensitive nature the dependence 
the energy gap across the (III) face, the choice the potential. 
Further, depends markedly the method employed. Howarth, when 
using the cellular method obtained energy gap 0.24 and when using 
the method augmented plane waves obtained energy gap approxi- 
mately 2.0 ev. reasonable conclude therefore that the accuracy with 
which may calculated these methods not great. 


-18.0 
awe 
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The position the peak the density states obtained Jones (1937) 
the center the zone face and obtained from optical data. the density 
states were calculated one the more elegant methods, then would 
appear that the accuracy with which the position the peak might deter- 
mined would directly related the accuracy with which might 
calculated. view the lack agreement between the methods and the 
difficulty concerning the choice potential, was felt that was worth while 
calculate the density states for silver using the method based the 
approximation nearly free electrons. 

The energy gap for silver was determined method suggested Jones 
(1934). the accuracy given the perturbation method, 

(4.1) AE@™ (2—f,)™ a*® 
where the atomic number and the atomic scattering power for rays 
sin and are the lattice parameters copper and silver 
respectively; the distance from the origin the (III) zone face. The 
scattering powers may obtained from tables (Peiser, Rooksby, and Wilson 
1955). The value for for copper chosen Jones from the optical data 
4.1 ev. For silver given (4.1) 4.8 which agrees moderately 
well with absorption edge silver beginning 4.0 (Mott and Jones 
1936). The density states calculation was carried out using value 4.8 ev. 
the lower value were used the peak would shifted 0.4 higher, but the 
form the curve would not appreciably altered. 

The results the calculation are shown Figs. and with the correspond- 
ing results for copper obtained Jones. noted that according 
this approximation the energy the Fermi surface occurs well beyond the 
peak for silver where the density states decreasing. the Fermi surface 
were spherical, the density states would not decrease until e/a 1.36 
(Hume-Rothery 1947). 

Since the explanation the experimental results presented this paper 
depends the density states decreasing with increasing electron atom 
ratio both silver and copper, appropriate consider briefly some the 
other experimental evidence bearing this point. Rayne (1957, 1958) 
measured the electronic specific heat and 
a-solid solutions. The results his measurements indicate initial rapid 
increase the density states with electron atom ratio followed 
levelling off even slight decrease. Jones (1957), however, has shown that 
when the Fermi surface lies close zone boundary the 
interaction leads term linear with temperature. suggested that this 
factor responsible for the apparent rapid rise the electron specific heat 
the brasses. Childs and Penfold (1957) criticizing the proposed explanation 
the observed variation the magnetic susceptibility with concentration 
the copper solid solutions (paper cited specific heat data indication 
that the density states increasing the Fermi surface. order explain 
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Fic. Density states for copper (Jones 1937) and silver function energy. 
Fic. Density states for copper and silver function e/a. 


the high ratio copper, silver, and gold, Klemens (1954) 
suggested that the Fermi surface touches the zone boundary all three 
metals. concluded from the work Smit (1952) and Mortlock (1953), 
who investigated the influence shear strain the thermoelectric power 
gold, and silver, and copper, respectively, that the Fermi surface touches the 
zone boundary all three metals. the other hand, Jones (1957) calculated 
expression for the thermoelectric power monovalent metals which, 
positive agreement with experiment, requires that large and 
positive, may just before the Fermi surface touches the zone boundary. 
Recently Rowland (1958) has measured the Knight shift the solid solutions 
formed silver with the Vth period elements. decrease the Knight shift 
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observed with increasing electron atom ratio. This may explained 
either decreasing density states increase the p-character the 
wave functions the Fermi surface. Cohen and Heine (1958) have recently 
surveyed the conflicting data this subject. Their attempt discredit the 
Jones (1937) model is, however, marred misunderstanding the model. 


DISCUSSION RESULTS 


paper was suggested that the atomic susceptibility the alloys 


where the mole fraction; the susceptibility the solvent metal 
ion; the susceptibility the solute ion with charge equal its valency 
plus the susceptibility the localized electrons where v—1. 
xr(x, the net paramagnetic susceptibility the conduction electrons, 
which function the electron atom ratio and the lattice parameter 
the solid solution a’. was assumed that 


where and are the density states the Fermi surface for lattice 
with lattice parameter and the lattice parameter the alloy respectively. 
The constant given 


(5.3) 


where the net paramagnetic contribution the conduction electrons 
the pure solvent, which equal the difference between the susceptibility 
the pure metal and its singly charged ion; the density states 
the Fermi surface the pure metal. Beyond the peak (see Figs. and 
the density states varies approximately linearly with e/a. The density 
states this region thus given 


The initial slope the susceptibility, concentration curve thus given 


where the proportional change the lattice parameter for unit concentra- 
tion (Owen 1947). 

order apply the method interpreting the results suggested 
first necessary decide value for Since the self-consistent field 
with exchange for has not been done and since the approximate method 
correcting for exchange used may unreliable this atomic number, 
alternative methods arriving value are necessary. which 
turns out small, neglected (5.5), then both and 
may calculated from the experimental slopes the Ag(IV) 
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series. assumed that xg+ approximately independent the solvent, 
then simultaneous equations may set the form 

where [dx the experimentally determined slope; 
and xzn+ and are the susceptibilities the singly charged 
and Gat solute ions which were calculated II. There are six sets equations 
this type. The results the six determinations and the probable errors 
Xagt and are and respectively. All 
susceptibilities are expressed c.g.s. e.m.u. per mole. 

Pauling (1927) calculated for semi-empirical method. 
Pauling’s empirical values agree within 10% with values calculated from 
self-consistent fields with exchange, where comparisons may made. The 
diamagnetic susceptibility the singly charged silver ion calculated from 
the self-consistent field (Black 1935) this corrected 
approximately for the effect exchange multiplying 0.8 then 
Xagt These two theoretical values for are reasonably 
good agreement with the experimental value. They not agree, however, 
with the value approximately obtained from inorganic salts 
(Pacault 1946). Pauling’s value for xag+ used the remainder this 
discussion. 

comparison may made this point between the theoretical value 
and the experimental value determined from the Ag(IV) series. The 
ratio the derivative the density states with respect the electron 
atom ratio, the density states found from the results 
Hence apparent good agreement with the experimentally 
determined value. must admitted, however, that the degree agreement 
may undetermined extent fortuitous. 

order demonstrate clearly the influence the solvents, Figs. 
plotted for each series solutes the two solvents. 
The results for the Cu(IV) systems have been taken from 
noted that the resulting curves fall into two series, the copper series and the 
silver series. referring (5.5), seen that the ordinate the sum 
the contribution the solute ion and the rate change the contribution 
from the conduction electrons. assumed that approximately 
independent the solvent, then the difference the ordinates given va- 
lency must due the difference the rate change the contribution 
the conduction electrons the two solvents. 

Since the self-consistent fields for and Sb+ have not been 
determined, the initial slopes the susceptibility, concentration curves may 
not calculated for the Cu(V) and Ag(V) series was done for the Cu(IV) 
series II. The difference the slopes the susceptibility, concentration 
curves for the same solute the two solvents may, however, calculated 


(5.6) 


| | 
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Initial slopes the susceptibility, concentration curves adjusted for the influence 
the solvent ion. 

Fic. Initial slopes the susceptibility, concentration curves adjusted for the influence 
the solvent ion. 
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using (5.5), subject the assumption that approximately inde- 
pendent the solvent. The calculated difference for the four series will then 
compared with the experimental difference and the degree agreement 
used determine the validity the initial assumptions. The difference 
slopes given 


order calculate the right-hand side (5.7) the following data has 
been used: (Pauling 1927); (Hartree 
and Hartree 1936); x&,(0), —0.17 
electron per atom per electron volt per electron atom ratio; 0.31 
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electron per atom per electron volt; —0.21; 0.23 (Fig. 5). 
paper has been assumed that slightly beyond the peak 
the density states. 

The observed and calculated differences the initial slopes are shown 
Table IV. The sum the first two terms (5.7), 
the next term and the last two terms taken together have 
maximum absolute value seen that the term involving 
the density states makes major contribution the calculated slope. 
(Better agreement, indicated the figures parentheses may achieved 
for the silver group value used.) The degree agree- 
ment between the observed and calculated values suggests that the form 
(5.7) correct for the polyvalent solutes and that the solute ion contribution 
fairly independent the solvent. important emphasize that the 
form (5.7) appears correct even though the details the calculation 

The treatment does not explain, however, why, for example, the absolute 
value initial slope the system greater than that the 
system, whereas the reverse true the corresponding 
copper base alloys. Also the results for the copper-silver solid solutions 
not fit the scheme suggested for the polyvalent solutes. Thus, for example, 
(5.5) would predict for the slope the system whereas 
the observed slope The experimental slope was corroborated 
temperature which the equilibrium concentration exceeded the solute 
concentration. 

expressed c.g.s. e.m.u. per mole and electrons per 
atom per electron volt, then for copper determined from (5.3) approxi- 
mately and for silver approximately 100. Although not possible here 
attach any special significance the relative magnitudes, interesting 
consider the origin the constants. 

believed that there are three terms which contribute the suscep- 
tibility the electron gas; the Pauli susceptibility xs, the 
susceptibility, (Mott and Jones 1936), and the interband contribution, 
(Adams 1953). The Pauli susceptibility given 
(5.8) 
where Bohr magneton. attempt was made show that the 
ratio the diamagnetism, the paramagnetism, —1/3 for the 
special form the relationship the approximation nearly free 
electrons. The integration the formula was done, however, 
over only one quarter the crystal, with the field the (111) direction, 
about which the function symmetrical. the direction the field and the 
(001) direction are made coincide and the integration done over the whole 


surface, then given 


where the free electron energy and the energy gap across 


HENRY AND ROGERS: MAGNETIC SUSCEPTIBILITIES 923 


the zone boundary. The contribution therefore, according 
the approximation nearly free electrons, equal constant times the 
density states. The sign positive, when decreases the 
system becomes less paramagnetic. Klee and Witte (1953) attempted 
explain the observed variation the magnetic susceptibility the ternary 
system using the paramagnetic term. Adams 
(1953, 1954) also pointed out that gives paramagnetic contribution 
when the Fermi surface not overlapped. Jones (1934) implied that 
always gives diamagnetic contribution, but this not true because, contrary 
his suggestion, not small the approximation nearly 
free electrons. 

Adams (1953) discussed the magnetic properties electron gas when 
the Fermi surface near zone boundary and found diamagnetic term, 
Which large under the same conditions and the same order 
magnitude. This contribution, proportional the density states 
(Adams 1959). Adams and Zitter (1954) interpreted the results Klee and 
Witte providing experimental demonstration the strength the 
interaction. They were able make their claim because, according 
their calculations, which responsible for the rapid rise the observed 
paramagnetism just before the Fermi surface touches the zone boundary, 
constant after the Fermi surface touches the zone boundary. Since varies 
directly with the density states, the constancy implies constancy 
the density states the same region. With the zone involved the 
work Klee and Witte rapid decrease, even more rapid than that found 
for copper and silver, the density states would expected for 
therefore difficult understand the constancy this region 
claimed Adams and Zitter. not fact constant, then becomes 
difficult distinguish between and 

From this discussion thus seen that may regarded the sum 
the constants proportionality relating xs, and the density 
states and not surprising that and should different. 

result the measurement the magnetic properties these alloys, 
something may said, general way, regarding the rate change 
the density states function electron ratio. The susceptibility, 
concentration curves show abrupt changes slope and some systems 
such and copper-gallium, the relationship quite linear 
electron atom ratio 1.30. reasonable expect that the 
passing through peak the density states curve would have noticeable 
effect the magnetic properties, follows that there peak the density 
states curve for either copper- silver-based alloys the range electron 
atom ratio investigated. 

There are two other physical properties which have been measured for 
these systems and with which some correlation might sought; these are 
the change lattice parameter produced at.% solute (Owen 1947) 
and the change resistivity produced at.% solute (Linde 
19326). relationship has been suggested between the volume changes and 


J 
> 
> 

’ 

’ 


924 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


the changes resistivity Blatt (1957). Owen found that the lattice dis- 
tortion was approximately proportional the valency the solute. The 
solvent effect the magnetic properties does not appear depend directly 
the relative rates change the lattice parameter with concentration. 
Thus found that the rate change the lattice parameter with concen- 
tration Ag(V) and Cu(V) Ag(IV) whereas the absolute 
values the adjusted initial slopes (Figs. and are Ag(V) Cu(IV) and 
Cu(V) There may small volume effect such suggested 
equation (5.2) but there evidence prove this. 

would like take this opportunity correcting error The 
atomic susceptibility the singly charged aluminum ion, was calculated 
Slater’s method and reported The correct value 
which agrees much better with which the 
sum and Shown Table II. 


SUMMARY 


The results the measurement room temperature the variation 
the magnetic susceptibility with concentration are presented for the solid 
solutions formed the four elements following silver the periodic table 
with both silver and copper solvents, and for the solid solutions, formed 
the four elements following copper the periodic table with silver 
solvent. These are compared with the results, obtained earlier, paper II, for 
the solid solutions formed the four elements following copper the 
periodic table with copper solvent. 

The four series are similar that the rate increase the diamagnetism 
the solid solutions series increases with valency. the silver series 
the rate greater than the copper series. both solvents the initial slope 
the susceptibility, concentration curves are approximately linear functions 
valency. 

The results the measurements are discussed terms model suggested 
paper The explanation appears give both qualitative and quantitative 
results reasonable agreement with experiment for the 
valent systems. 

The explanation suggested requires that the density states the Fermi 
surface both copper and silver decreases upon the addition polyvalent 
solute. less specific alternative this would require only that the 
expression for the initial slope the susceptibility, concentration curve con- 
tains term which proportional the excess valency, and that the 
constant proportionality depends only the solvent. This was done for 
the silver base alloys with the [Vth period solutes and the constant pro- 
portionality determined. 

order test the explanation offered paper II, was necessary 
calculate the density states for silver the approximation nearly 
free electrons (Jones 1937). found that the Fermi energy corresponding 
electron atom ratio 1.0 occurs energy about 1.0 electron 


ee = 
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volt beyond the peak. The calculated ratio the rate decrease the 
density states, beyond the peak, copper silver, 0.8. 

argued that since the susceptibility, concentration curves show 
abrupt changes slope and are, some cases, quite linear, there peak 
the density states curve either copper- silver-based alloys with 
electron atom ratio the range 1.0 approximately 1.30. 

The solid solutions the copper-silver alloy systems are anomalous and 
behave though they are phase mixtures, i.e. the susceptibility approxi- 
mately linear function the susceptibilities the pure metals. 
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ILLUSTRATION TYPE TRIPLE CORRELATION 
MEASUREMENT WHICH CAN EASILY 


ABSTRACT 


The reaction has been used illustrate the simplifications 
introduced the interpretation triple angular correlations choosing 
target and bombarding particles zero spin and observing the emitted 
particles, this case neutrons, counter fixed the beam. The angular 
correlations the gamma rays with respect the incident beam then depend 
only upon the properties the final states the residual nucleus. The angular 
correlation the electric quadrupole 2.03-Mev gamma ray predicted uniquely 
theory and this prediction has been verified experimentally. The angular 
correlations the 1.28-Mev and 2.43-Mev gamma rays have yielded for the 
amplitude mixing ratios +0.25+0.05 —3.4+0.5 and —0.26+0.08 

—1.10+0.16 respectively. addition, the experiment provides illustration 
the value the recently discov ered technique gamma-ray 
discrimination organic scintillator. 


INTRODUCTION 


This paper provides illustration particular type triple angular 
correlation measurement which can considerable value the deter- 
mination nuclear spectroscopic data from nuclear reactions. this type 
measurement aligned nuclear states are selected the detection the 
second radiation with counter which axially symmetrical with respect 
the incident beam direction. The expression for the angular correlation 
the third radiation, gamma ray for example, then the same that derived 
for the decay nucleus aligned low temperatures (Cox and Groot 
1953). the case the selection aligned nuclear states from nuclear 
reaction the relative contributions from the different magnetic substates are 
unknown unless the details the formation and decay the compound 
nucleus are specified. However, the magnetic substates that can populated 
are limited the spins the incident and emergent particle and the target 
(Litherland 1959). The example given this paper illustrates that 
certain cases the summation over the magnetic substates can eliminated. 
The gamma-ray angular correlations then depend only upon the properties 
the states the residual nucleus. 

the experiment described the reaction was chosen 
for study because the zero spin the target nucleus and the bombarding 
particle.* neutrons from the reaction are then observed the beam 
only the +1/2 substates the levels the residual nucleus are 
populated. This because the orbital angular momentum vector the 

received March 22, 1960. 
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*In cases where the spins the target nucleus and projectile are not zero, sometimes 


possible remove the uncertainties interpretation due the summation over the magnetic 
substates (Litherland and Ferguson 1959). 
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neutron perpendicular its direction motion. Consequently, only the 
spin the neutron can have projection the direction motion the 
neutron. The full theoretical expression for the angular correlation under 
these conditions given Section this paper. 

The reaction has number other advantages which 
make suitable for triple correlation measurements the type described 
above. 

(a) The states have been studied number reactions and the 
spins and parities the ground, the 1.28-, 2.03-, and 2.43-Mev states have 
been assigned (Endt and Braams 1957). addition, some measurements 
Litherland and Gove (1958) the gamma rays from the He? reaction 
have given information the cascading from the 3.07- and 3.62-Mev states 
The most prominent gamma rays observed from the reaction 
are shown Fig. and the yield curves these gamma 
rays are shown Fig. 


~15 Mev EXCITATION 


3.07 (32+ 52+) 
2.03 
1.28 
O+ ‘e+ 
26 


Fic. The low-lying states 


GAMMA-RAY ENERGIES | t } GAMMér RAY ENERGIES 
1.28 Mev 2.03 Mev 
1.79 Mev 1.59 Mev 


| 


4.0 5.0 4.0 5.0 


ALPHA~PARTICLE ENERGY IN Mev ALPHA-PARTICLE ENERGY IN Mev 


RELATIVE INTENSITY 


Fic. Yield curves the gamma rays from the reaction obtained 
Litherland and Gove (1958). The energy scale should increased approximately kev. 


The low Q-value the reaction +0.017 Mev (Endt 
and Braams 1957) ensures that low bombarding energies only limited 
number low-lying states are excited. Furthermore, the competing 
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reaction has the unfavorable Q-value —3.19 Mev which 
implies low yield for this reaction the He* energies used (i.e. between 
and Mev). Consequently, the gamma rays from the low-lying states 
can readily studied with little interference from other reactions 
from the interference multiple cascading from high excited states. 
the experiment described this paper, was possible measure both 
the angular correlation the gamma rays coincidence with the neutrons 
and also, some cases, measure the angular distribution the gamma 
rays with respect the He* beam with the neutrons unobserved. 

(c) The low Q-value the reaction implies that the neutrons leaving 
excited states will fairly low energy Mev). This feature 
the reaction relevant the correction for the finite size and asymmetry 
the neutron counter used since the correction smaller for low angular 
momentum neutrons than for high angular momentum neutrons. This effect 
discussed briefly Section this paper. 

The practicability the experiment described this paper depends largely 
the neutron gamma-ray discrimination properties some organic scin- 
tillators (Brooks 1959; Owen.1959) because necesssary use efficient 
fast-neutron detector which insensitive gamma rays, for the coincidence 
measurements. method discriminating against gamma-ray pulses from 
large liquid scintillator described Sections and Sections and 
this paper are devoted discussion the experimental results and 
the conclusions and possible future applications the type triple correlation 
measurement discussed this paper. 


THEORETICAL CONSIDERATIONS 


Theoretical expressions for the angular distribution the gamma rays 
with respect the incident beam alpha particles, suitable for the analysis 
the reaction have been given Sharp al. (1954). 
Expressions for the angular correlation the gamma rays observed coin- 
cidence with the neutrons are also given. the latter case, the expression 
given quite complex and involves extensive summations over vector coupling 
coefficients. shown Litherland and Ferguson (1959) considerable simpli- 
fications arise the neutrons are observed axially symmetrical counter 
180° the incident beam. this case, the theoretical expression 
for the observed angular correlation the gamma rays from the reaction 


Q,P; (cos 6). 


equation the angular correlation the gamma rays ex- 
pressed terms Legendre polynomials. P(m) the population the 
mth substate the nuclear state spin which decays another state 
spin emitting gamma-ray quantum total angular momentum 
The coefficient has been tabulated Sears 


. 


930 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


and Radtke (1954). The coefficient has been tabulated 
Sharp al. (1957). The attenuation coefficients can calculated for 
cylindrical gamma-ray counter the method suggested Rose (1953). 
Tables these coefficients for 5in. diameter long, diameter 
6in. long, and diameter long Nal crystals have been pre- 
pared Gove and Rutledge (1958) and Rutledge (1959). The full expression 
for the angular correlation must include sum over the different multipole 
radiations that may present. This expression, for the 
case, then given equation 


equation the amplitude the electric quadrupole radiation 
present divided the amplitude the magnetic dipole radiation. ex- 
pression similar equation together with table angular distribution 
coefficients has been given Wapstra al. (1959). 

Equation also applies the case when the neutron unobserved and 
the angular distribution the gamma rays measured with respect the 
beam. this case, the populations the gamma-ray emitting state have 
regarded unknown parameters determined experimentally. If, 
however, the angular correlation the gamma rays are observed coinci- 
dence with the neutrons detected infinitely small counter situated 
the incident beam then only the +1/2 substates are populated 
discussed earlier. The angular correlation the gamma rays from the 
reaction then depends only upon the properties the states the residual 
nucleus. This conclusion does not depend upon the assumption that the 
reaction proceeds through isolated resonance the compound nucleus. 
For example, the angular correlation electric quadrupole gamma ray 
arising from 5/2 1/2 transition uniquely predicted 


The 2.03-Mev gamma ray which was studied the experiment reported 
this paper arises from such transition. Since the spins the other states 
studied are known, measurement the angular correlation coefficients 


oz | 

=~ | 
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VALUE 
=. GAMMA RAY 
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GAMMA RAY 


MIXING RATIO 


Fic. theoretical values the triple correlation coefficient are shown 
function the amplitude ratio, together with the experimental values (discussed 
Section obtained from the angular correlations the 1.28- and 2.43-Mev gamma rays. 


LITHERLAND AND McCALLUM: TRIPLE CORRELATION MEASUREMENT 931 


the de-excitation gamma rays serves determine the magnetic dipole 
electric quadrupole mixing ratio The angular correlation coefficient 
shown function Fig. for the case 3/2 1/2 transition. 

Since infinitely small neutron counters cannot used detect the neutrons 
the beam, necessary consider the effect the finite size 
neutron counter which large enough give significant coincidence counting 
rate. this effect has been treated some detail Litherland and Ferguson 
(1959), only their conclusions will quoted here. the particular reaction 
discussed this paper the effect using small axially symmetrical neutron 
counter half angle introduce contribution the angular corre- 
lation the gamma rays from the +3/2 substates. the case the 
angular correlation the 2.03-Mev gamma ray discussed this paper the 
angular correlation then becomes 


equation and are unknown constants which are related the 
formation and decay matrix elements the compound nucleus. These con- 
stants can calculated only assumptions are made regarding the formation 
and decay the compound nucleus. For example, consider the case 
resonance formed the capture d-wave alpha particles 
this resonance then decays the emission neutron the 5/2+ 
2.03-Mev excited state can shown, considering the triple 
angular correlation detail, that the coefficients and are +9/4 and 
—3/8 respectively. Figure shows the exact calculation the angular 
correlation coefficients function the angle (no longer small). Two 
cases are shown; the compound state spin decays either the emission 
neutron can result only the population the +1/2 substates 
the residual nucleus there attenuation the correlation coefficient 
result the finite size the neutron counter. 


APPARATUS 


Experiments were carried out using mya 4.5 5.0 Mev ions 
from the Chalk River 3-Mv electrostatic accelerator. The ions bombarded 
thin metallic target (50 kev thick for 4.7-Mev alpha particles) deposited 
onto 0.02-in. tantalum backing.* The gamma rays were detected 
in. diameter long Nal crystal, the front face which was 
6.5 inches from the target. The crystal could rotated from 120° 
with respect the beam. The neutrons were detected 


*Obtained from the Electromagnetic Separator Group, Atomic Energy Research Establish- 
ment, Harwell, England. 
from Nuclear Enterprises Ltd., Winnipeg, Manitoba, Canada. 
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(a) 
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The effect the finite size the neutron counter illustrated particular case. 


discrimination circuit the type first described Owen (1959) was used 
discriminate between neutrons and gamma rays. 


The experimental arrangement illustrated Fig. the beam 


the liquid scintillator could not placed position ensure complete 
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axial symmetry because the scintillator did not completely fill its container. 
was therefore placed the beam axis with its axis symmetry perpendicular 
the beam distance inches from the target. The effect the deviation 
from axial symmetry the reaction studied small for reasons discussed 
below. 

block diagram the electronic equipment used the experiment 
shown Fig. The Nal crystal was optically coupled Dumont 
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Fic. Block diagram the electronic equipment used the experiment. 


LINEAR 
AMPLIFIER 


CHANNEL 
KICKSORTER 


6364 photomultiplier and Dumont 6363 photomultiplier was coupled 
the liquid scintillator. Two pulses were taken from the multiplier connected 
the Nal crystal, pulse for pulse height analysis from the 
last but one dynode and pulse for coincidence analysis from the 
anode. Three pulses were taken from the multiplier connected the liquid 
scintillator, pulse from the fifth dynode from the anode, 
pulse from the third dynode from the anode, and space-charge limited pulse 
from the last dynode. The space-charge limited pulse can used dis- 
criminate between neutron and gamma-ray induced pulses from the liquid 
scintillator the manner described Owen (1959) and Litherland 
al. (1959). The Dumont 6363 multiplier was used with approximately 1500 
volts between anode and photocathode. The last dynode-to-anode voltage 
was made remotely variable between the limits approximately and 
volts. For the measurements described this paper and the particular 
photomultiplier used the optimum last dynode-to-anode voltage was found 
volts. Under these conditions, the neutron-induced pulses the last 
dynode (in the energy range studied) had prominent negative component 
whereas the gamma-ray induced pulses were predominantly positive. After 
amplification, neutron-induced pulses could selected discriminator 
which was triggered the negative portion the space-charge limited pulse. 

*The pulse height, from proton-induced pulse liquid scintillator, not course 
linear function the proton energy. 
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The outputs the neutron discriminator and fast coincidence unit 
the type described Bell a/. (1953) (resolving time musec) were 
applied slow coincidence unit which consequently gave output only 
when neutron and gamma ray were fast coincidence. This output pulse 
was used open linear gate which transmitted 100-channel transistor 
pulse-height analyzer* those gamma-ray pulses which were coincidence 
with neutrons. second pulse height analyzer displayed the spectrum all 
neutrons within energy range determined the settings single- 
channel pulse-height analyzer and which were coincidence with output 
from the discriminator set the space-charge limited pulses. The scaler 
shown Fig. was used monitor that the coincidence spectra 
each angle could normalized given neutron flux. 


EXPERIMENTAL PROCEDURE 


When were setting the electronic equipment shown Fig. the 
experimental procedure adopted was first adjust the neutron gamma-ray 
neutrons and gamma rays. This was most readily accomplished observing 
the space-charge limited pulses oscilloscope which was triggered the 
pulses from single channel analyzer set the linear pulses. The single 
channel analyzer window was adjusted with the help calibrating gamma-ray 
sources such cover the approximate range proton recoil pulses 
interest the experiment. this way was possible adjust the last 
voltage that the neutron and gamma-ray induced pulses 
were opposite sign over the proton recoil energy range studied. Adjusting 
the level the discriminator the space-charge saturated pulses and throw- 
ing the switch shown Fig. alternately the positions marked 
and gamma and permitted the discriminator level 
which gamma-ray pulses were eliminated rapidly achieved. During the 
experiment, the operation the gamma-ray discriminator was 
continuously monitored the oscilloscope. 

The satisfactory operation the complete electronic system was ascertained 
observing neutron gamma-ray coincidences from the source and 
then observing the lack coincidences from source. Finally, small 
adjustments were made using the neutrons and gamma rays from the 
reaction. 

For convenience setting the apparatus, the upper and lower limits 
the proton recoil spectrum were defined the discriminator set the 
space-charge limited pulse. This possible because the negative part the 
space-charge saturated pulse, characteristic neutron-induced pulses, first 
increases and then decreases with increasing proton recoil energy. Conse- 
quently, the neutron-energy discriminator bias was adjusted that neutrons 
feeding the fourth excited state (shown Fig. were not counted. 
This setting was checked the absence 1.79-Mev gamma rays from the 
decay the fourth excited state the first excited state The action 


*Obtained from Computing Devices Canada Limited, Ottawa, Ontario, Canada. 
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the gamma-ray discrimination circuit illustrated Fig. 
where the spectrum gamma-ray pulses from the Nal crystal, fast 
coincidence with pulses from the liquid scintillator, shown (a) with the 


100 


1.79 Mev ~2.03 Mev 


50 


100 + 


COINCIDENCE COUNTS 


GAMMA-RAY ENERGY Mev 


Fic. Coincidence pulse spectra taken without neutron gamma-ray discrimination; 
and with neutron gamma-ray discrimination. 


discriminator inoperative and (6) with the discriminator operation. The 
1.79-Mev peak seen curve (a) arises from the de-excitation the 3.07 level 
This peak not visible curve (6) showing that (i) the neutron 
gamma-ray discriminator bias set sufficiently high eliminate co- 
incidences from the excitation the 3.07 level and (ii) coincidences are 
effectively discriminated against. This latter feature necessary because the 
3.07-Mev state also cascades through the 2.03-Mev state emitting 
1.04-Mev and 2.03-Mev gamma rays. further check the rejection 
coincidences the gamma-ray discrimination circuit 
coincidence spectrum was measured both with and without the discrimi- 
nation circuit operative. reduction coincidence counting rate factor 
over 200 was obtained with the discrimination circuit operation. 

The angular correlations were carried out 4.74- and 4.94-Mev bom- 
barding energy with the neutron counter fixed first and then 90° 
the beam. The energy the ions was measured relative the 4.00-Mev 
resonance the reaction (Gibbons and Macklin 1959). The 
angular correlations the 1.28-, 2.03-, and 2.43-Mev gamma rays were 
obtained simultaneously 4.94-Mev bombarding energy. 4.74 Mev the 
yield 2.43-Mev gamma rays was too low permit accurate angular 
distribution made. 

Since the gamma rays traverse the neutron counter when they are observed 
near correction was made for the resulting attenuation. This correction 
was determined observing the angular distribution the 2.185-Mev gamma 
rays from radioactive source placed the position the target. The 
neutron counter was for this measurement. the correction was large 
(amounting ~40% between and 20°) would advantage detect 
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the outgoing neutrons near 180° and thus obviate the correction for attenuation 
the gamma rays. This possibility discussed below. 


RESULTS AND DISCUSSION 


Fig. the angular correlations the 2.03-Mev electric quadrupole 
gamma rays from the decay the 5/2+ second excited state the 


THEORY 


THEORY 


400 
300 
ro) 
© 200 


100 
Mev 
°° 30° 60° 90° °° 30° 60° 90° 
ANGLE @ ANGLE @ 


Fic. The triple angular correlations from the reaction leading 
excitation the level are shown for two incident alpha-particle energies and 
for two geometrical arrangements the counters. 

(a) The neutron counter situated the beam (Fig. and the coincidence counting 
rate observed function the angle the gamma-ray counter. 

(b) similar correlation with the neutron counter situated 90° the beam. 


1/2+ ground states are shown. The theoretical correlation, derived from 
equation for the case with infinitely small neutron counter fixed 
(marked (a) Fig. and infinitely small gamma-ray counter angle 
are corrected only for the finite aperture the gamma-ray counter. The 
calculations Gove and Rutledge (1958) were used for this correction. Agree- 
ment with theory quite good despite the finite size the neutron counter 
and its lack exact axial symmetry. The results least squares analysis 
the angular correlation data the two energies are given Table The 


TABLE 


Angular correlation coefficients 


energy 
4.74 Mev +0.62+0.06 
4.94 Mev +0.55+0.06 


Theory 0.53 


attenuated theoretical coefficients are also quoted Table The least 
squares fits reported this paper were carried out using the Chalk River 


Datatron computer. 


100 
™ 


AS 


LITHERLAND AND McCALLUM: TRIPLE CORRELATION MEASUREMENT 937 


The good agreement with theory, without the correction for the finite size 
and partial asymmetry the neutron counter, can explained partly the 
smallness the correction term given Section this paper and partly 
assuming that the neutrons the two resonances are predominantly 
s-wave, because shown Fig. for s-waves the predicted triple correlation 
independent the angle the neutron counter. This assumption supported 
the measurement, 4.74-Mev bombarding energy, the angular distribu- 
tion the 2.03-Mev radiation with the neutrons unobserved. These angular 
distribution coefficients are consistent with the neutrons being 60% s-wave 
and 40% d-wave. However, the observed triple correlations with the neutron 
counter 90° are different the two resonances studied and are different 
from the correlations showing that there are differing admixtures the 
higher angular momentum neutrons the two energies. 

Using the theoretical expression for the angular distributions given 
equation together with the measured angular distribution coefficients for 
the 2.03-Mev radiation given Table Ila, possible determine the 


TABLE 
Angular distribution coefficients 


energy Gamma-ray energy 
(a) 4.74 Mev 2.03 Mev —0.16+0.07 
(b) 4.94 Mev 2.43 Mev —0.48+0.08 


populations the substates the 2.03-Mev (5/2+) state These are 
quoted Table These population numbers give support the con- 


TABLE III 


P(+1/2) P(+3/2) 
(a) 4.74 Mev 2.03 Mev 0.02+0.03 
(b) 4.94 Mev 2.43 Mev 0.09+0.04 


jecture made earlier that the neutrons observed each resonance are pre- 
dominantly low angular momentum. 

The arguments outlined above are given because the experimental conditions 
were not ideal. Ideally, the neutron counter should have been axially sym- 
metrical and should have been adjustable distance from the target. The 
effect the finite size the counter could then have been investigated 
varying the distance the counter from the target and observing the change 
the angular correlation pattern. Such experiments would require axially 
symmetrical counter with hole through allow the beam reach the 
target. The distance the neutron counter would not then limited the 
necessity for the gamma-ray counter observe the gamma rays near 0°. Such 
experiments are being planned. order interpret the results discussed 
this paper, will assumed that the effects the finite size the neutron 


4 
Population numbers deduced from angular distribution coefficients 
r 
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counter are negligible compared with the errors determining the experi- 
mental coefficients. The measurements the angular distributions the 
2.03- and 2.43-Mev gamma rays and the triple angular correlation the 
2.03-Mev gamma ray can used support this assumption. 

The correlations the 1.28-Mev gamma rays obtained 4.74- and 4.94- 
Mev bombarding energies are shown Fig. and least squares fit the 


4.94 Mev 


= = 4.74 Mev 1 
o rT 
200+ ¢ ¢ 


T 


L 
60° go" 


ANGLE 6 


Fic. Triple angular correlations the 1.28-Mev gamma ray are shown for two incident 
alpha-particle energies. The neutron detector was the incident beam both cases. 


data expression the form (cos was obtained. From the 
two correlations mean value —0.05+0.07 was obtained. Fig. 
the experimental value the (cos coefficient given together with the 
theoretical value plotted function the magnetic dipole electric 
quadrupole amplitude mixture. The theoretical coefficients were obtained 
from equation with P(m) 1/2 for and P(m) for 
The theoretical value the coefficient has been modified for the effects the 
finite solid angle subtended the gamma-ray counter. The mixing ratio 

4.74 Mev the yield the 2.43-Mev gamma rays was too low permit 
accurate triple angular correlation measurement. The correlation the 
2.43-Mev radiation 4.94 Mev, with the neutron counter 0°, shown 
Fig. together with the correlation obtained with the neutron detector 
90°. least squares fit the data, with the neutron counter 0°, yielded 
value —0.75+0.07 for the (cos coefficient. The angular distribution 
the gamma ray with neutrons unobserved was measured three 
angles from which value for the (cos coefficient was 
deduced. This coefficient included Table and from the populations 
listed Table were deduced. the case the formation the 2.03- 
state these populations imply that the neutrons the 2.43-Mev state 
are low angular momentum. Consequently correcting only for the angular 
aperture the gamma-ray counter, mixing ratio —0.26+0.08 
was obtained from Fig. 


CONCLUSIONS 


Triple angular correlations from the reaction have been 
used determine values for the states The ambiguity the results 


| 
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4.94 Mev 


600— 


COINCIDENCE COUNTS 


Fic. correlations the 2.43-Mev gamma ray with neutron detector (a), 
and 90° for bombarding energy 4.94 Mev. The solid curve least squares fit 
the experimental data taken with the neutron detector 0°. 


can principle eliminated observing the linear polarization the 
gamma rays. This was done McGowan and Stelson (1958) for number 
gamma-ray transitions produced Coulomb excitation. addition, the 
experimental results described this paper have illustrated the considerable 
advantages gained the interpretation triple correlation data 
choosing zero-spin target and projectiles and observing the neutrons 
the beam projectiles. The simplification the interpretation is, however, 
obtained the price the experimental difficulty using counters which 
are axially symmetrical with respect the beam and which are situated 
180° the beam. Though such ideal counters were not used the 
experiment described above, the effect the finite size the counter and its 
lack axial symmetry was quite small because the neutrons observed were 
mainly low angular momentum. Since, general, these simplifying condi- 
tions not hold, further experiments are planned using neutron counter 
which axially symmetrical and which has hole through its center that 
neutrons can observed 180° the incident beam. This counter will 
remove the uncertainties due both the attenuation the gamma rays 
traversing the neutron detector when mounted the incident beam, 
and the azimuthal asymmetry the detector. addition, would also 
permit the experimental evaluation the correction the data for the 
finite size the neutron counter. 

Since and induced reactions can excite large variety nuclear 
states, attempt has been planned use axially symmetrical junction 
counter with hole through the center detect the alpha particles and 
protons from these reactions 180° the incident beam. junction counter 
with these characteristics has been successfully made and tested McKenzie 


and Bromley (1959). 
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MAGNETO-OPTICAL ROTATION VERY HIGH 


RICHARD STEVENSON 


ABSTRACT 


Magneto-optical rotation transmission through reflection from solids 
examined the classical free electron theory, with the view taking such 
measurement using fields the megagauss range. general the rotation 
markedly non-linear function the magnetic field, and some cases can change 
sign the field increases. For very low fields the rotation varies directly with 
but the high field limit the rotation varies inversely with the field. For 
substances which the intercollision time the electron small, measurements 
the Kerr rotation (i.e. reflection) will give the electron mobility function 
the magnetic field, and thus will give important data which can used 
conjunction with high field magnetoresistance experiments. 


INTRODUCTION 


Measurements the Faraday rotation microwave infrared 
frequencies have been suggested and used number times obtain values 
the effective electron mass non-ferromagnetic substances. The theory 
for low fields infrared frequencies outlined Mitchell (1955) indicates 
that the rotation has dependency and varies directly with the field. The 
theory such experiments microwave frequencies has been considered 
Rau and Caspari (1955) and more recently Lin (1959). More general use 
Faraday rotation experiments has been considered Stephen and Lidard 
(1959). Measurements Suhl and Pearson (1953), Stern and Myers (1958), 
and Brown and Lax (1959) indicate that acceptable technique for 
measuring the effective electron mass low fields. 

When the magnetic field very large, approaching the megagauss mark, 
appears that the energy density the field can affect the band structure 
the solid. Since some cyclotron resonance experiments are done high fields, 
would seem worth while have another type experiment which could 
performed under similar conditions and would give similar information about 
the effective electron mass. Then would possible study what effects, 
any, could attributed the presence the very strong field itself. For 
example, the high field experiments Keyes al. (1956) show variation 
effective mass, the reason for which somewhat obscure without subsidiary 
information. 

OUTLINE THE THEORY AND RESULTS 
The solid described optical terms complex index refraction 
and the index for radiation circularly polarized one sense distinguished 

received August 31, 1959. 
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from that the opposite sense the extinction coefficient not zero 
the light transmitted the Faraday effect reflected the Kerr effect will 
elliptically polarized and the major axis the ellipse rotated through 
angle terms the indices and N_, the rotations and ellipticities 
(Argyres 1955) are: 

(a) Faraday rotation and ellipticity 


(1) 


Kerr rotation and ellipticity 


N,N_-1 
where and refer the ‘‘real and imaginary parts the frequency 
incident radiation, the thickness the sample, and the velocity 


light. 
The rotation the polarized light due the electrons the solid, which 
have the equation motion 


where displacement vector the electron, 
mean collision time the electron, 
the resonant frequency bound electron, 
incident radiation field, 
external magnetic field. 


All quantum mechanical effects are contained attributing the electron 
m*. The term involving will dropped, that the 
theory strictly applies only free electrons. The type analysis has been 
explained Sommerfeld (1954) and its applicability has been considered 
Givens (1958). 

The incident beam plane-polarized radiation represented by: 


(4) 
Then using equations can obtain 


where have defined the static conductivity, 


N,—N_ 


ill 
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with being the number electrons per unit volume, and the cyclotron 
frequency 


The relation between the a-c. conductivity and the static conductivity 
derived somewhat similar manner, Dresselhaus a/. (1955). Accord- 
ing our definition (4) the incident radiation, obtain 


And from these equations the terms (1) and (2) can evaluated. 


(7) 
where 


Then from (1) formula can obtained for the Faraday rotation, which gives 
the general dependence the magnetic field 


where the electron mobility 


has been introduced. When the external field small, i.e. (8) reduces 
either the formulae derived Mitchell (1955) Rau and Caspari (1955) 
depending whether not large small. 

The very high field behavior some interest. Even low fields the rota- 
tion becomes non-linear with respect has been reported Brown and 
Lax (1959). the field increases the Faraday rotation reaches maximum, 
and can change sign. For extremely high fields the rotation loses its depen- 
dence the effective electron mass and tends vary inversely with the field 


Thus, however useful this technique may low fields obtain values 


n 
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m*, does not appear that such experiment can used with the very high 
pulsed fields give information complementary that obtained from 
cyclotron resonance. 

Measurements the Kerr rotation evidently will give useful information 
about electron transport metals high fields. For most metals the inter- 
collision time the electron small, the order sec less ordinary 
temperatures. Then experiment most metals will satisfy the conditions 
that infrared frequencies and fields megagauss. From the 
factors (7) one can obtain expression for the Kerr rotation metals 


eT 
m* 


which formula contains the assumption that the static conductivity the 
metal large. give idea the size the effect can consider metal 
such zinc with the representative parameters sec and .02 
which easily measurable. 

The interesting feature (10) that allows measurements the electron 
mobility function the magnetic field made independently 
magnetoresistance measurements. Thus knowing the effective electron mass, 
possible obtain data the field dependence the electron collision 
time This topic some theoretical interest which there are evidently 
experimental data. The measurement electron mobility should also 
useful the accumulation data high field behavior transport properties. 

Since one effective mass has been used the model strictly correct only for 
solids with spherical energy surfaces, which is, course, consistent with the 
application free electron theory. The quantum mechanical treatment 
this type problem gives the same results the free electron theory, and 
specific high field behavior which causes depression the Fermi level and 
removal the degeneracy the electrons (Kahn and Frederikse 1959) does not 
invalidate the use here the free electron model. some substances may 
possible observe oscillatory behavior the optical properties the high 


fields. 
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DIFFUSION PHENOMENA AND ISOTOPE EFFECTS THE 
EXTRACTION FISSION-PRODUCT XENON AND 
KRYPTON FROM IRRADIATED 


ABSTRACT 


investigation the diffusion fission-product xenon and krypton from 
irradiated powder function temperature has revealed isotopic 
fractionation and also double-valued activation energy. The apparent fission 
yields and show abnormal enrichments factor 10. These 
enrichments appear related the precursor half-lives xenon. When 
isotopic fractionation exists, the diffusion results exhibit extremely low energy 
activation. 


INTRODUCTION 


the determination fission yields using mass spectrometer techniques, 
frequently desirable extract both the solid and gaseous fission products 
from the same sample. The fission gases may extracted either heating 
the irradiated samples high temperatures, dissolving them completely. 
During study the fission yields nuclides produced the thermal fission 
sample fission gas was obtained simply heating irradiated 
sample finely divided plutonium oxide under conditions such that only part 
the fission-product gas was recovered. The yields the xenon isotopes 
observed for this sample were found differ markedly from those obtained 
when the total fission gas was extracted dissolution. 

The fractionation the xenon isotopes which occurred the partial extrac- 
tion the fission gas was too large accounted for the basis simple 
diffusion theory applied xenon gas. Also, the isotope effects for the various 
xenon isotopes were not mass dependent. This rather puzzling result, coupled 
with the growing interest the use oxides potential fuel-element material, 
suggested more complete study should made this effect. 

Experiments were therefore designed determine the xenon and krypton 
isotopic abundances fission gas fractions extracted various temperatures 
from neutron-irradiated powders. From these experiments, the isotope 
effects could followed and the diffusion activation energies determined. The 
results are reported below. 


EXPERIMENTAL 


Sample Preparation 

prepare uranium oxide, analytical grade uranyl nitrate was decomposed 
heating air about 750° chemical analysis showed the prepared 
material After this heat treatment, the oxide was ground until 
would pass 250-mesh screen (30 micron particle radius). Samples 
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grams were then transferred into cylindrical quartz vials having inside 
diameter 5.5 mm, outside diameter 7.5 mm, and length mm. 
The quartz containers were evacuated, sealed under vacuum, and then placed 
into self-serve aluminum capsules for subsequent irradiation the N.R.X. 
Reactor Chalk River. Three samples were irradiated, one after the other, 
the same self-serve position, for periods 2.5, and days. The neutron 
flux this position was about neutrons 

portion the uranium oxide was separated means sedimentary 
tube order determine the particle size distribution. This measurement 
indicated that the most probable radius was the range from 
microns. 


Rare Gas Extraction 

Following the irradiation, the samples were stored for months permit 
the decay the major activities. the end this period, the quartz vials 
were placed into high-vacuum lines similar that shown Fig. The system 


CALCIUM 


FURNACE 


SAMPLE 
TUBE 


DIFFUSION 
PUMP 


TRAP 


FURNACE 


Fic. schematic representation the system used extract the fission-product rare 
gases various temperatures. 


was thoroughly evacuated the vial was opened. With the vial opened, 
the system was connected the vacuum pumps through charcoal trap held 
liquid air temperature. Such arrangement ensures collection any 
fission-product krypton and xenon released from the sample. furnace which 
had been thermally calibrated surrounded the sample region. The sample was 
raised through series temperatures, remaining each value for hours. 
The rare gases liberated each temperature plateau were collected and 
transferred calcium furnace for purification. The purified gas was then 
collected into sample tube and removed from the line for subsequent mass 
spectrometric analysis. 


— 
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Mass Spectrometry 

The isotopic analysis the rare gas samples was performed with 10-in. 
directional-focussing mass spectrometer which employed electron multiplier 
for the detector. The resolving power the instrument was 650 which 
sufficiently high permit the separation the rare gas isotopes from their 
respective isobaric hydrocarbon impurities. high sensitivities these hydro- 
carbons are always present within the mass spectrometer tube itself and would 
mask small rare gas samples lower resolution were used. The isotopic 
abundances normal krypton and xenon were found agree within with 
the accepted values given Nier (1950). The over-all sensitivity the 
instrument permits the analyses gas samples having volume 


RESULTS 


Effects 

The isotopic abundances and the corresponding volumes the fission- 
product rare gases were accurately determined for each temperature plateau 
and estimate was made the corresponding volume released gas. The 
irradiation conditions and the xenon abundances for the three samples 
are given Table Ten spectrograms were obtained for each sample every 


TABLE 
The abundance fission-product xenon evolved from function temperature 


Isotopic abundance 


Tempera- Sample Sample Sample 
ture irradiated 2.5 days irradiated days irradiated days 
plateau, 


— 


Note: The error for all abundances is less than 2%. The recorded values have been normalized to Xe1%*, with 
a value of 1.00 for the Xe! yield. The recorded values all exceed this value because of the contribution from the 
reaction Xe!35 (n,y) Xe136, 


temperature plateau and resulted most probable error less than for 
the relative abundances. each case, the observed abundances were nor- 
malized should noted that the yield changes for each 
sample because the contribution function the irradia- 
tion time. The abundances observed for samples extracted between 800° and 
1200° were found excellent agreement with the fission yields 
(Blades 1956) for these isotopes. These samples contained the bulk 
the fission gas and fractionation effects were not detectable. Normal xenon, 
which can detected the presence was observed for the 30° 


5 
0.29 0.49 0.76 1.02 0.39 0.60 
150 1.70 0.76 0.80 1.02 2.50 0.83 
300 2.69 1.02 0.87 1.02 2.55 1.41 
500 0.65 0.84 0.99 1.02 0.76 0.90 1.12 
800 0.38 0.59 1.00 1.02 0.38 0.58 1.12 
1200 0.36 0.55 1.00 1.02 0.36 0.55 1.12 
Expected 
Ss 
n 
S 
a 
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fractions but dropped rapidly zero the temperature was increased. 
This indicates that small amount normal xenon was adsorbed the 
surface the particles during the preparation the samples. 

The data tabulated Table show effects which would not expected 
priori. The and abundances not remain constant relative 
but show remarkably high enrichments the gas fractions extracted 
150, 300, and 500° For example, the abundance for sample enriched 
plateau. With the exception sample the abundance remains 
constant within the experimental error for all temperatures. These effects are 
most clearly seen the graphical presentation given Fig. where the 


SAMPLE | SAMPLE 2 SAMPLE 3 


RELATIVE FISSION YIELD 


% 300 600 900 12000 30 600 900 12000 300 600 900 i200 


TEMPERATURE 


Fic. Histograms showing the observed fission yields function extraction tempera- 
ture for the three samples studied. The true expected fission yield always observed 
1200° The dashed curves indicate the Maxwellian distribution which best fits the enrichment 
factor. 


observed abundance plotted function temperature. should noted 
that the abundances observed the region from 800° 1200° are identical 
with the accepted values for indicating detectable isotope fractionation. 

The experimental results yield only histogram and, therefore, con- 
venient find smooth curve which gives the differential abundance every 
temperature within the interval. inspection the histograms for and 
shown Fig. suggests that the enrichment factor can most easily 
represented Maxwellian distribution. The Maxwellian parameters were 
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obtained determining the area and first moment the respective tempera- 
ture-abundance histograms. The curves determined this manner are shown 
the dashed lines Fig. The fit very good and suggests that even higher 
enrichments could obtained measuring the abundance fission gases 
released within smaller temperature intervals. 

The krypton isotopes released from the samples different tempera- 
tures showed deviation from their normal fission pattern. Because krypton 
more abundant nature and less abundant fission than xenon, was 
impossible make meaningful measurements the fission-product contribu- 
tion the low temperature fractions. Reliable measurements isotopic 
abundance could only made from 500° 1200° where the normal 
krypton contamination small compared with the fission-product krypton. 
Therefore, can concluded that 500° too high temperature observe 
any appreciable krypton fractionation, if, indeed, such effect does exist 
for the krypton isotopes. 

The cumulative volume gas released and including each temperature 
plateau recorded Table These values have been normalized the total 


TABLE 
The cumulative volume rare gas evolved for 3-hour heating periods 


Temperature, 


150 300 500 800 1200 


Total xenon 0.13 1.00* 
Total krypton 0.12 


*All the fission-product rare gas is assumed to escape when the U;0Os is maintained at 1200° C 
for 3 hours. The values recorded represent an average of the three samples investigated. The 
total volume of fission gas produced in sample 2, is calculated to be 6.3 X 1074 cc for xenon and 
1X 10-5 ce for krypton (both at N.T.P.). 


volume gas which assumed have escaped the time the powder 
has been maintained temperature 1200° for 3-hour period. This 
assumption appears valid since the measured volume 1200° was good 
agreement with the volume calculated from the known sample size and irradia- 
tion conditions. The values quoted represent the average the three samples 
for all isotopes. These data show that the abnormal xenon enrichments involve 
only about 0.5% the total volume and also that from 500° 1200° the 
xenon and krypton diffuse the same rate. 

The results given Tables and permit one determine the cumula- 
tive abundances, i.e., the abundances which would observed all the gas 
released temperature was collected and examined. Such calculation 
was carried out and the results are given Table III. interesting note 
abundances. 

obtain more complete information about the enrichment mechanism, 
subsidiary experiment was conducted. For this experiment, samples 
were prepared manner described previously; however, only particles 
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TABLE III 


The cumulative abundance fission-product xenon evolved from function 
temperature 


Isotopic abundance 


Tempera- 
ture Sample Sample Sample 

0.29 0.49 0.76 1.02 0.39 0.60 0.97 1.08 0.53 0.71 1.00 1.12 

150 0.40 0.51 0.76 1.02 0.51 0.69 0.99 1.08 1.29 0.98 1.12 

300 0.84 0.61 0.78 1.02 1.63 1.09 1.01 1.08 1.70 1.09 1.01 1.12 

500 0.81 0.65 0.82 1.02 0.85 0.92 1.02 1.08 0.67 0.78 1.01 1.12 

800 0.53 0.61 0.94 1.02 0.41 0.61 1.01 1.08 0.38 0.58 1.01 1.12 

1.00 1.02 0.36 0.55 1.00 1.08 0.36 0.55 1.00 1.12 


1200 0.36 0.5 


Note: The error for all abundances is less than 3%. The recorded values have been normalized to Xe!%, 


having radii from microns were used. This particle-size selection was 
achieved the use sedimentary tube. The selected was then uni- 
formly distributed within powdered-sulphur matrix order provide 
trap for any fission fragments recoiling from the uranium particles. Two such 
samples were then irradiated low neutron flux avoid melting the sulphur. 
one sample, the rare gases trapped within the sulphur were removed 
subliming the sulphur from the mixture under vacuum. The solid fission 
products the second sample were obtained treating the mixture with 
system and extracting from the water phase. both cases, 
least 10% the fission products were found the sulphur matrix. 
The isotopic abundances the fission-product krypton, xenon, cesium, and 
neodymium from the sulphur matrix were determined. The krypton and 
neodymium did not deviate more than 10% from their fission yields; 
however, the xenon and cesium yields, which are given Table IV, showed 


TABLE 


Mass number 


131 132 133 134 135 136 


137 


Xenon ratio 1.0 
Cesium ratio 1.7 


*For a ‘“‘normal case"’ the recorded value would be 1.0. 
+Cs!*3 and Xe!*4 have been normalized to 1.0. 


marked enrichments. The results presented Table are the ratios 
(observed yield)/(true yield) with renormalization the and 
ratios equal 1.0. With such representation the data, value 1.0 will 
obtained enrichment present. Since cesium monoisotopic (viz., 
only lower limit for the normal cesium contamination can given. 
Therefore, the and ratios represent only lower limits since they 
are measured relative which has been assumed entirely radio- 
genic. The sulphur experiment reveals xenon fractionation not unlike that 


The relative quantities* xenon and cesium appearing sulphur matrix 
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for and, addition, discloses similar effect for cesium. The data also 
indicate that major isotope enrichment occurs the krypton and neo- 
dymium mass regions. 


Diffusion Effects 

The properties the gross diffusion process, whereby the rare gases are 
released from the particles, can ascertained rather simple considera- 
tions. Since the particles are approaching sphericity, the problem can 
reduced one diffusion from sphere. obtain solution the problem, 
the diffusion equation 


must solved with the boundary conditions, 


(i) uniform initial concentration time 
(ii) zero surface concentration for 


The solution equation (1) under these conditions 


n=] 


where the particle radius and the initial concentration. The quantity 


Defining F,, the fraction the total gas released time can 
write 


where the particle volume. For small values i.e., equation 
(4) reduces 


solving either equation (4) (5) the appropriate regions, can 
obtained function the parameter Table gives experimental 
values, therefore, means this calculation, the corresponding Dto/R? 
value can determined. The results this computation are given Table 
where the parameter given for each temperature plateau studied. 

the diffusion process can described equation the form 


then plotting against the inverse temperature straight 
line should result, the slope which the activation energy graph 


*The term can written since are constant throughout the experiment. 


| 

(9) R’ 
S 
33 


952 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


TABLE 


The diffusion parameter function temperature for the 
fission-product rare gas and system 


Xenon Krypton 
Temperature, 

150 

800 0.13 0.12 
1200 1.00 1.00 


*The quantity Ft represents the fraction of the total gas released after a 
time fo at a given temperature 7. For the calculations in this table it has been 
assumed that had the sample been heated initially to temperature T, all the gas 
released at lower temperatures would be evolved at temperature 7. Presumably 
the correct value to use would lie between the integral and differential quantities 
evolved. However, since the two limiting cases make no significant change in the 
results, the former assumption has been retained. 


the data plotted this manner shown Fig. immediately apparent 
from Fig. that two activation energies exist and also that the process can 
satisfactorily described equations the form given equation (6). 


2000 600 600 500 400 300 

Fic. relationship found for the diffusion constant and The dimensionless 
quantity proportional since remained constant throughout the experiment. 
existence two activation energies indicated the different slopes found above and 


below 600° 
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temperatures below 410° the predominant term has activation energy 
3000 cal (0.13 ev) while that above 410° has value 39,000 cal 
(1.7 ev). This effect appears hold only for xenon since the krypton 
activation energy remains 39,000 cal over the region studied. The 
curve Fig. can described the equation 


Since sec and lies between and microns, the value 
the frequency factors and can calculated. When this done, the 
errors the exponent represent the spread the radius the particles. 
interesting that the fractionation was observed for gas samples extracted 
below 500° which precisely where the low activation process prevalent. 


DISCUSSION RESULTS 


The most striking feature this investigation has been the large isotope 
effect found for Xe™! and and the association this with extremely 
low activation energy. rather early propose mechanism for this 
effect; however, several interesting features are revealed the data. 

connection with isotopic effects for rare gases, have noted that there 
little fractionation apparent among the krypton isotopes for 
and These isotopes, which exhibit little fractionation, all 
share the common property having short-lived precursors. marked 
contrast this, and have relatively long-lived precursors 


This apparent correlation fractionation with precursor half-life suggests the 
diffusion properties and concentration distributions tellurium and iodine 
are quite different from those xenon. Another possibility that the fractiona- 
tion related the recoil range the fission fragments. This seems rather 
unlikely, the basis our knowledge recoil ranges (Finkle al. 1951). 
Also, for the case powder, the initial distribution the fission 
products expected isotropic regardless the range except for the 
boundary between the powder and the quartz container. 

consider only the correlation fractionation and half-life, then would 
appear that there tendency for tellurium and iodine concentrate the 
surface the particles greater extent than for xenon. When the 
iodine decays, there enrichment Xe™! and the surface 
grain boundaries. One then expects large enrichment and low activation 
energy associated with these isotopes. 

The results the sulphur matrix experiment depends upon recoil ranges 
and surface concentrations. restrict ourselves the latter, then interface 
diffusion tellurium and iodine and the recoil xenon following decay 
may possibly account for the isotopic distribution found the sulphur matrix. 
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The cesium yields associated with the sulphur fraction also exhibit isotopic 
effects. Here one finds that the cesium enrichment increases with decreasing 
half-life the xenon precursor. This suggests that cesium, iodine, and tellurium 
all behave similar manner the surface the particles and all 
behave quite differently from xenon. The xenon appears have relatively 
low mobility and does not concentrate the surface grain boundaries 
the particles. 

Although the diffusion data obtained the course this work are qualita- 
tive nature, they do, however, suggest the following interesting features. 
The apparent existence two activation energies can associated with the 
observed isotopic effects. The activation energy below 600° depends upon 
the particular isotope studied, while above this temperature, all the xenon 
and krypton isotopes have approximately the same value. Thus, below 600° 
the apparent activation energy could related the precursor half-life 
the xenon. The existence this double value may possibly explain the wide 
range activation energies found the literature (Stubbs al. 1957; 
Lindner and Matzke 1959). 

The results obtained this work have suggested many avenues investiga- 
tion which should explored complete understanding the observed 
effects achieved. With this mind, active program has been initiated 
and results should shortly forthcoming. 
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DIATOMIC 


KovAcs 


ABSTRACT 


Explicit expressions are obtained for the intensity distribution the branches 
and bands. They are valid with satisfactory accuracy for 
any value the coupling constant the state. 


Theoretical investigations intensities electronic bands have hitherto 
mostly dealt with transitions which the orbital angular momentum remains 
unchanged changes +1. sure, some authors have extended their 
investigations also other cases but these are very special kinds (Schlapp 
1932; Van Vleck 1934). Recent experimental data (Herzberg 1959), however, 
make necessary the study further cases particularly transitions between 
and states the same and different multiplicity. the present paper 
the and transitions are considered. 

well known that the interaction neglected, the elec- 
tronic part the wave function any state the molecule can separated 
into the product two factors, one depending the orbital motions, and 
the other the electron spins. account interaction the wave 
functions are modified and, result, the above-mentioned transitions are 
permitted with small intensity. 

Let denote the electronic parts the wave functions the individual 
components the A... states, neglecting the interaction 
taking into account the interaction, the perturbed wave functions 
will linear combinations the unperturbed wave functions the same 
and multiplicities From the properties the operators and 
readily verified that the interaction joins terms for which 
+1, +1, and One the perturbed wave functions the 


following: 


where constant order unity, and are constants much less 
than unity if, shall suppose, the interaction small. 

Taking account the degeneracy the states which differ only the 
orientation the resultant relative the axis the molecule (+2 and 
and completing the wave functions the parts arising from the rotations 
and vibrations the nuclei, the new wave functions will the following: 

Contribution from the Department Atomic Physics, Polytechnical University, Budapest, 
Hungary. 

the explicit form these wave functions, see, for example, Kovacs (1958). 

tThese have form similar the factors inside the parentheses eq. (5). 
Can, J. Phys. Vol. 38 (1960) 
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il 


where and are the vibrational and rotational wave functions the 
molecule. Since are concerned with definite vibrational level each 
the two states question, the vibrational part the wave function merely 
introduces into the intensities factor independent and will therefore 
omitted what follows. The task theory the calculation the 
so-called factors for all branches occurring the rotational transitions. For 
this the expressions for the matrix elements the components the electric 
moment the molecule 


where the angle een the and preferential direction space, 
taken the z-axis, may write the following explicit expressions for the 
amplitudes 


where the exponent takes the values case the and 


whether are dealing with states. The constants are inde- 
pendent for example given the following expression: 
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While these constants determine the possibility transitions between the 
above-mentioned states, need not consider them for the problem the 
intensity distribution except case the transition where two 
constants independent each other occur, that only one these may 
eliminated. The absolute values the amplitudes are 


sin 


cos 


(6) 


They may found paper Kronig (1927). 

The squares (4a) and (4c) multiplied and summed over the magnetic 
quantum numbers give the factors referring the transitions 
and (a) and these are found the third column Tables 
and III. 

However, multiplet states always belong case Therefore 
for the state introduce wave functions the form 


3 


(7) 
= J+1 
Hence the amplitudes become 
(8) 


obtained the basis (8) are listed the third column Table II. 

The states can general represented well the formulae Hund’s 
case (a) only for small values the rotational quantum numbers. With 
increasing rotational numbers there transition toward case 
The wave functions valid with satisfactory accuracy for any value the 
coupling constant can then represented 


1 
where the transformation matrix elements are the following: 
= 
(10) 4 IC. So, IC. (J) 


(J) = VCi( J) 


Sig-1 


| | 
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Here for normal states positive) 


while for inverted states negative) the coefficient 
instead 

obtain the factors the intermediate case between case (a) 
and for the term, the amplitudes must written the following 
form: 


where From (12) the total intensities the transitions 
the fourth column Tables II, and III. Table III the constant 
given (see (4c)). The expressions are valid for any value the 
coupling constant and contain special cases the expressions for the tran- 
sitions 


X X - 3A inv (@), Xx - 3A(d) 


when The expressions for the latter special case are given 
the fifth column Tables II, and the tables the terms were 
the state but and therefore for the transitions 
branches should replaced the designation (for the first 
index) and the Tables and III necessary interchange the branches 

Comparing the intensities each the branches the transition 
with those the corresponding branches the transition found 
that for low rotational quantum numbers the intensities the branches 
those the branches are greater whilst for high rotational 
numbers the intensities the branches are smaller, those the 
are greater than the case the transition having the same coupling 


*The exact expression for the following Ci(J) 
satisfactory. similar formula holds the case 

(1932) has dealt with the transition but since had interchanged two 
columns the transformation matrix the term, his results are not correct. 

3, 
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constant This comparison helps recognize the kind transition since 
the method the missing lines the neighborhood the zero gap not 
always applicable because too much overlapping. 


The author deeply indebted Dr. Herzberg for suggesting the problem. 
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NOTES 


THE RELATIVE ABUNDANCES COSMIC RAY NUCLEI CHARGE 


The charge spectrum heavy nuclei the primary cosmic radiation has 
been the subject large number investigations. recent years stacks 
nuclear emulsions have been exposed very high altitudes, and result 
general agreement has been reached that the flux light nuclei 
not negligible the top the atmosphere. However, disagreement still 
exists the exact abundances the various groups nuclei, and much 
this uncertainty due the extrapolation the experimentally determined 
fluxes the top the atmosphere. 

For the present results, stack Ilford stripped emulsions, in. 
um, was flown from Brownwood, Texas (geomagnetic latitude, 
March 1958, for 7.5 hours ceiling height 136,000 feet. 
The residual pressure for this height 2.7 With the packing material 
surrounding the stack the total amount matter above the emulsions was 
equal 3.2 g/cm?, which less than previous investigations. Less extra- 
polation therefore necessary obtain the fluxes the top the atmosphere. 

All the emulsions were scanned along line below the top edge for 
tracks with grain densities least six times that minimum ionization. The 
tracks selected for measurement had have minimum projected length 
per emulsion sheet and zenith angle 80°. order estimate the 
efficiency scanning for tracks about one quarter the stack was scanned 
experienced observer. comparison with the relative numbers tracks 
found other scanners makes believe that tracks were missed. 
distinguish the tracks from those slowing down a-particles all the 
suspected tracks were followed until they either interacted left the stack. 
total 612 tracks satisfying the selection criteria were found. 

Since the cutoff energy geomagnetic latitude 41° about 1.5 
Bev/nucleon, all the primary nuclei entering the stack are relativistic. 
principle one ionization measurement should therefore sufficient for 

charge determination. the present results, however, both the 6-ray density 
and gap density each track with nuclear charge were determined. 
least 150 four grains more and 600 gaps specified depth 
regions the emulsion were counted each track. These measurements were 
carried out two observers working independently, one using Koristka 
and the other Leitz microscope. The results each observer were cali- 
brated with the aid charge-indicating interactions, and was found that the 
calibrations were exactly the same for both observers, whereas the gap 
*Issued N.R.C. No. 5687. 


tNational Research Council Postdoctorate Fellow. Now National Physical Research 
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densities differed about 10%. This was due the different optical con- 
ditions the two microscopes. The results each observer were, however, 


self-consistent. 

Figure shows the combined normalized results the two observers. The 
charge resolution good Nevertheless, several tracks falling 
near the borderline between the boron and carbon groups have been cross- 


CHARGE FROM GAP COUNTING 


Fic. Comparison the charge determination density and gap density measure- 
ments. The crosses the figure represent tracks for which the measurements could not 


carried out the specified depth regions the emulsion. 


checked both observers. significant differences were found. Table 
gives the numbers light 5), medium 9), and heavy 
10) nuclei different zenith angle intervals. They are also expressed 
percentages the total numbers nuclei each interval. Table gives the 
same results corrected for loss tracks due interactions occurring the 
emulsions above the scan line and also for particles which entered the stack 
during ascent and descent (Waddington 1956). Finally, Table III gives the 
total numbers light, medium, and heavy nuclei the interval 60° 
corrected the top the atmosphere using the Bristol fragmentation probabi- 
lities published Rajopadhye and Waddington (1958). These numbers are 
also given absolute fluxes and percentages. 
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TABLE III 
Fluxes nuclei different charge groups the top the atmosphere 


Light nuclei Medium nuclei Heavy nuclei 


Total number nuclei (0° 60°) 111.5 351.9 108.3 


Absolute flux particles/m? sec sterad 1.14+0.11 3.60+0.19 
Relative abundance (%) 20+2 19+2 


The value 0.20+0.02 obtained for the relative abundance light nuclei 
experimental agreement with the results Bristol (Waddington 1957), 
Chicago (Koshiba 1958), Rochester (Engler a/. 1958), Minnesota 
(Freier al. 1959), Washington (O’Dell 1959), and Torino (Garelli 
1960). the other hand, inconsistent with the value about 0.05 
obtained Appa Rao al. (1958) and the value 0.30 Noon 
(1957). The absolute fluxes are considerably lower than those obtained 
previous observers stacks flown over Texas. This could perhaps 
attributed the fact that the present flight took place during period 
low cosmic ray intensity. 

are now following all the primary and secondary tracks until they 
either interact leave the stack. Full details the present work, including 
the fragmentation probabilities obtained, will published later date. 


Appa Rao, K., Biswas, S., DANIEL, R., NEELAKANTAN, A., and 1958. 
Phys. Rev. 110, 751. 
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TRIBOLUMINESCENCE MERCURY BUBBLER 


Weak light emission may observed when various gases are allowed 
bubble through depth several millimeters mercury contained pyrex 
glass quartz vessel atmospheric pressure. The glowing region appears 
consist numerous small flashes light occurring several centimeters 
above the surface the liquid depending the value the gas flow rate. 
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Using glass transmission grating spectrograph and f/4 quartz 
prism instrument, the spectra associated with the observed luminosity has 
been photographed and identified. 

The production luminosity did not appear depend critically the 
design dimensions the mercury bubbler. The two types bubbler used 
are shown Fig. The brightest emission occurred when helium and argon 


Fic. forms mercury bubbler. 


were employed the bubbling gases. the experiments with helium, the 
most intense spectral features were found due 0.02% air impurity 
initially present the gas.* These spectra consisted the bands 
relatively high intensity. Weaker features due transitions from upper levels 
above were also identified. The bands were missing when argon was 
used the bubblant, otherwise the emission corresponded that observed 
for helium. 

Wavelengths and relative intensities these spectra are given Table 
Intensities are visual estimates plate blackening arbitrary scale 
600. intensity indicates faint trace after 5-hour exposure with 
the spectrograph used with 103 film. intensity 600 associ- 
ated with the strongest feature which can just photographed minutes 
with the f/4 spectrograph employing 103 plates. The spectral features 
photographed with the quartz optics are illustrated Figs. and 2c. With 
nitrogen the bubbling gas, the second positive system well 
spectra and appeared; however, the emission was weaker than for 
either helium argon. With ethane the bubblant, the resulting luminosity 
was too feeble for spectra obtained conveniently. 

should like thank Dr. John Herron N.B.S. for the mass spectrometer analysis. 


(a) (b) 
| ! 


PLATE 


0,0 


(c) 


4000 3000 2500 


Mercury bubbler and Tesla discharge spectra. 

(a) Mercury bubbler with helium; exposure hours; f/4 quartz spectrograph. 

(b) Tesla discharge through helium atmospheric pressure; exposure minute, minutes, 
minutes; quartz spectrograph. 

Mercury bubbler with argon; exposure hours; quartz spectrograph. 
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NOTES 


TABLE 
Spectra observed the mercury bubbler with helium bubblant 


Identification 


Intensity (v’, Other 


2953 

2962 

2977 

3089 
3136 
3159 
3285 
3309 
3339 
3371 
3500 
3537 
3577 

3672 

3710 

3755 

3805 

3943 

3999 

4047 Hg (8S; 3Po) 

4059 

4200 

4270 

4344 

5016 

5770 

5790 

6678 


NoTE: m = masked. 
“Tabulated wavelengths (Pearse and Gaydon 1950; Harrison 1939), 


The observed emission example triboluminescence initiated the 
hydrodynamic interaction between flowing gas and liquid. the mercury 
bubbler, electrostatic charging mercury droplets known take place 
means series complex processes inherent bubble formation and 
disintegration (Harper 1957). the present case, the luminosity almost 
entirely confined the region near the vessel wall (most clearly seen the 
bubbler Fig. although pressures few millimeters mercury where 
the breakdown potential the gas reduced, the glowing region seen 
spread throughout the volume the vessel manner characteristic 
electrical discharge process. appears that discharge between charged 
droplet and the surrounding silica surface the source the luminosity 
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described above. electrical discharge between the liquid and the vessel 
wall result frictional charging the liquid due its mechanical motion 
relative the wall also thought significant process. The hypothesis 
discharge mechanism the source luminosity further supported 
Tesla coil excitation the same mixture atmospheric pressure. 
Except for the occurrence the NO, system, the spectra this source (Fig. 
26) correspond the spectra obtained from the mercury bubbler. Since colli- 
sional deactivation processes are likely important atmospheric pressure, 
the observed decrease luminescence intensity with increasing complexity 
the flowing gas probably reflection the relative deactivation efficiencies 
associated with each the gaseous species employed. 

The observation bands the bubbler with helium the bubblant 
but not with argon implies that metastable rare gas atoms act inter- 
mediate species the transfer energy the nitrogen molecule. order 
excite radiation from the state least 18.8 are 
15.6 for the ionization and 3.2 for the excitation the resulting 
ion. The required energy not available from collisions involving metastable 
readily supplied collision with helium the long lived level lying 
above the ground state the atom. 


Harrison, 1939. wav elength tables (John Wiley Sons, Inc., New York). 

(Chapman Hall, Ltd., London). 


RECEIVED 1960. 
NATIONAL BUREAU STANDARDS, 
D.C. 


LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 


previous the first day the month issue. proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
reference standard (with assumed frequency 192 631 770 c.p.s.). Frequency devia- 
tions from nominal are quoted parts per 10". negative sign indicates that the frequency 


below nominal. 


GBR, 16 ke/s 


Date, 

April MSF, 54-hour 24-hour WWVB 
1960 average* average ke/s 

3 — 162 —172 — 167 N.M. 

5 — 164 — 167 — 166 —-79 

Y — 143 — 167 — 169 —81 

8 —173 — 164 — 162 -—79 

9 — 163 — 166 — 163 N.M. 
167 —160 —82 
14 * —170 — 166 — 162 N.M. 
N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M. N.M. N.M. N.M. 
N.M N.M. N.M. N.M. 
N.M N.M. N.M. 
N.M N.M. N.M. 
28 —141 — 148 —144 N.M. 
Midmonthly mean — 159 — 163 — 160 -79 

Midmonthly mean 
WWV 


Note: N.M. no measurement. 
*Time of observations: 00.30 to 04. 
11.00 tol 
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Experimentally Measured Vibrational Levels 


electron selector has been used obtain experimental values for the separation 
some vibrational levels NO*, and These are believed the first such measure- 
ments made using electron bombardment techniques. The results are particular 
interest because spectroscopic data exists for this, the simplest molecular structure. 

The lower curve Fig. represents the potential energy diagram the neutral hydrogen 
molecule. Using Morse function for the potential, the levels the various vibrational 
states may calculated quantum mechanical methods (e.g. Herzberg 1950). These 
theoretical values have been verified spectroscopic observations (e.g. Richardson 1934). 
The waves the various levels represent qualitatively the probability the nuclear separa- 
tion. ordinary mass spectrometer source temperatures, the population the first excited 
level neutral hydrogen negligible, the molecules being the ground state 

The upper curve Fig. represents the potential energy diagram This has been 
calculated theoretically for Morse function (Hylleraas 1931) and the vibrational levels 
established. These values agree with those calculated extrapolation the observed 
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Fic. 1. Vibrational energy states of He and H*. 
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vibrational levels for various excited states the neutral hydrogen molecule (e.g. Richardson 
1929). However, direct experimental measurements have been made these levels. The 
values nuclear separation are indicated qualitatively for several levels wavy 

Appearance potential measurements the ion formed from electron bombard- 
ment have been made using the electrostatic electron selector (Marmet and Kerwin 1960). 
Several typical curves are shown Fig. where the ion current has been recorded 
function the electron energy. Since the energy spread the selector beam about 0.03 ev, 
relatively sharp breaks are obtained, and number them are seen separated 
about 0.25 ev. 


Fic. 2. Appearance potential curves for H}. 


Average values the differences electron energy required produce successive breaks 
were taken from runs. They are here compared with the differences between the first six 
vibrational levels calculated from Richardson’s constants (Herzberg 1950). 


SELECTOR THEORETICAL 
v=1 0.269 ev 
A 0.27+0.02 ev 2, 0.253 ev 
B €.25+0.02 ev 3, 0.237 ev 
C 0.23+9.02 ev 4, 0.221 ev 
D 0.21+0.02 ev 5, 0.208 ev 


6, 0.192ev 


Our selector values clearly show the narrowing vibrational intervals. While our experi- 
mental error does not permit unambiguous assignment, the intervals appear correspond 
the first vibrational levels. 

may noted Fig. the most probable transition from the ground state the mole- 
cule the ion not the first vibrational level, but one such that probable nuclear 
separation equal that the molecule principle). There certain 
probability for transition nearby levels. Experimentally, the probabilities may calcu- 
lated from the slopes the appearance potential curve segments. our case, the relative 
probabilities were: 


This confirms the general picture Fig. and the fact that the nuclear separation greater 
for the ion than the molecule. our first level indeed the ground one, then the most probable 
transition lies 0.75 above the theoretical 15.4-ev jump from the molecular ground the 
ionic ground state. Such value 16.1 has been frequently observed, and our values 
agree with Bleakney’s analysis the conflicting data the appearance potential 
(Bleakney 1932). 
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placed. Manuscripts and illustrations should carefully checked before they are submitted. 
Authors will charged for unnecessary deviations from the usual format and for changes 
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and the principal findings, required, except Notes. 
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numbered, and typed after the text. The form the citations should that used current 
issues this journal; references papers periodicals, titles should not given and 
only initial page numbers are required. The names periodicals should abbreviated 
the form given the most recent List Periodicals Abstracted Chemical Abstracts. All 
citations should checked with the original articles and each one referred the text 
the authors’ names and the year. 

Tables.—Tables should numbered roman numerals and each table referred the 
text. Titles should always given but should brief; column headings should brief and 
descriptive matter the tables confined minimum. Vertical rules should not used. 
Numerous small tables should avoided. 
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are appear the reproduction should ruled black ink. Paper ruled green, yellow, 
red should not used. All lines must sufficient thickness reproduce well. Decimal 
points, periods, and stippled dots must solid black circles large enough reduced 
necessary. Letters and numerals should neatly made, preferably with stencil (do NOT use 
typewriting) and such size that the smallest lettering will not less than high 
when the figure reduced suitable size. Many drawings are made too large; originals 
should not more than times the size the desired reproduction. Whenever possible 
two more drawings should grouped reduce the number cuts required. such groups 
drawings, large drawings, full use the space available should made; the ratio 
height width should conform that journal page in.), but allowance must 
made for the captions. The original drawings and one set clear copies (e.g. small 
photographs) are submitted. 

Photographs.—Prints should made glossy paper, with strong contrasts. They 
should trimmed that essential features only are shown and mounted carefully, with 
rubber cement, white cardboard, with space between those arranged groups. mount- 
ing, full use the space available should made. Photographs are submitted 
duplicate; they are reproduced groups one set should mounted, the duplicate 
set unmounted. 
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